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ERRATA 

1.  Page  6,  3rd  paragri^h,  line  7.  Delete  "The  'correction'"  and  sub¬ 
stitute  "Ilf  an  observation  station  reports  a  height  and  a  wind,  the  correction" 

2.  Page  6,  3rd  paragraph,  line  8.  Delete  "grldpolnt"  and  substitute 
"location  of  the  observation  (Interpolated  from  the  four  surrounding  grldpolnts)" 

3.  Page  6,  3rd  paragrjyih,  line  9.  Delete  "winds,  and"  and  substitute 
"winds.  This" 


4.  Page  6,  3rd  paragr{4>h,  line  9.  Change  "the  equation"  to  read  "the 
corresponding  equation" 


5.  Page  6,  Eq.  (2-1).  Cha.nge  to  read 

mkfv. 


C  =  W(d) 


1  +  n 


g 


mkfu^ 


Z  . 
x,yj 

(2- la) 


6.  Page  7,  5th  line  from  bottom.  Delete  "The"  and  substitute  "Ilf  a 
station  reports  a  height  only,  the  correction  equation  Is 
rS(x  ,y  )  +  nZ  .  I 

c  =  W(d)  [  (2-lb) 

Ilf  a  station  reports  a  wind  only,  Eq.  (2- la)  Is  used,  but  the  terms  nZ^j^g 
are  omitted.  7  The" 


7.*  Pages  17  and  21.  Substitute  revised  Figs.  4-l(d)  and  4-5  provided. 


^Result  of  machine  error. 
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8. *  Pages  79  and  80.  Substitute  the  following  equations  for  45— 55‘N: 

250  mb  Zq  =  588.77  +  0.96044  (300Z^)  +  1.7419  (300Tjj) 

250  mb  Tq  =  -  511-33  +  1.0568  (250Z^)  -  1.0568  (300Zjj)  -  0.96757  (300Tjj) 

200  mb  Zq  =  734.04  +  0.95107  (250Zp)  +  1.9341  (250T^) 

200  mb  Tjj  =  -  464.99  +  0.78515  (200Z^j)  -  0.78515  (250Zjj)  -  0.78136  (250Tjj) 

9. ')'  Page  105.  Substitute  revised  portion  of  table  provided. 

’"Result  of  machine  error. 


vlii 


#• 


The  revised  material  on  this  and  the  following  page 

may  be  cut  out  and  pasted  over  the  original  printing.  ^  4!  ®  "S 

O  ^  ^  CM 


t  Revised  Fig.  4-1  (d). 

Cut  out  and  paste  onto  page  17. 
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Revised  portion  of  table.  ■  » 
Cut  out  and  paste  onto  page  105. 
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Revised  Fig.  4-5  is  on  the  following  page. 
Cut  out  and  paste  onto  page  21. 
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ABSTRACT 


^Vertical-extrapolation  regression  equations  for  use  in  specifying  initial- 
guess  fields  for  objective  analyses  at  100  mb  were  derived  (for  four  midseason 
months)  through  use  of  a  screening  procedure.  Tests  of  the  equations  on 
independent  data  proved  the  equations  to  be  stable,  and  the  results  for  mid¬ 
season  months  are  considered  excellent.  Further  tests  of  the  equations  on  data 
from  one  month  on  each  side  of  each  midseason  month  indicated  that  approxi¬ 
mately  95%  of  the  equations  are  useful  and  stable  for  the  "side"  months. 

An  e.xperiment  comparing  the  regression  equations  derived  here  with 
similar  equations  developed  by  the  Navy^pi^^in  the  200— 100-mb  extrapolation 
interval  illustrates  that  the  newly  derived  equations  yield  significantly  lower 
root-mean-square  errors  than  the  Navy  equations, 

A  build-up  analysis  procedure  incorporating  the  vertical  extrapolation 
equations  is  described,  and  the  results  of  limited  development  testing  of  the 
analysis  procedure  are  discussed. 

Results  of  an  experiment  that  compared  analyses  employing  initial  guesses 
given  by  12-hr  persistence  with  analyses  employing  the  build-up  procedure 
indicated  that  the  build-up  produces  significantly  lower  root-mean-square  errors. 
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1.0  INTRODUCTION 


Objective  analysis  of  tenq)erature8,  heights,  and  winds  at  stratospheric 
constant-pressure  surfaces  is  a  complex  problem  because  of  relatively  sparse 
data  coverage  and  increased  instrumental  and  radiation  errors  in  the  available 
data  [1,  12].  Obviously,  a  good  initial  guess  is  essential  to  any  numerical 
objective  analysis  at  stratospheric  levels  because  over  the  no-data  and  sparse- 
data  areas  the  "analysis"  will  effectively  be  the  initial  guess.  Persistence  has 
been  suggested  [6]  as  an  initial  guess,  but  it  is  not  adequate  over  oceans  and 
other  sparse-data  regions  in  which  previous  analyses  are  also  questionable 
because  of  lack  of  data.  The  major  aim  here,  therefore,  has  been  directed 
toward  generating  the  best  possible  initial-guess  fields  for  height,  temperature, 
and  winds  for  100-,  50- ,  and  30-mb  objective  analyses  over  the  hemisphere. 

Recently,  Teweles  and  Snidero  discussed  some  of  the  problems  involved 
in  objectively  analyzing  stratospheric  constant-pressure  surfaces.  Their 
paper  (13)  describes  a  statistical  build-up  of  a  pattern  from  a  200-mb  analysis 
to  100  mb,  using  vertical-extrapolation  regression  equations  developed  by  the 
Navy  [4] .  The  Navy  equations  use  only  the  height  and  temperature  at  200  mb 
for  extrj^olation  to  100  mb.  Regression  equations  were  developed  by  the  Navy 
for  extrapolation  of  height  and  temperature  through  four  other  intervals; 

200  to  50  mb,  100  to  50  mb,  100  to  30  mb,  and  50  to  30  mb.  In  each  case,  the 
height  and  temperature  of  the  lower  level  were  the  only  predictors  found  to  be 
useful.  Several  parameters  were  tested  as  possible  secondary  predictors. 

These  parameters  were  various  combinations  of  squares  and  products  of 
temperature  and  height  at  the  base  of.  the  extrapolation  interval  and  the  next 
lower  level,  and  the  difference  between  the  temperature  at  the  extrapolation 
base  and  the  temperature  at  the  next  lower  level.  It  was  found  that  the  inclusion 
of  these  did  not  significantly  improve  the  basic  regression  equations. 

One  possible  reason  that  the  secondary  predictors  selected  did  not  improve 
the  results  over  the  basic  regression  equations  is  that  the  coefficients  for  the 
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primary  predictors  were  held  constant.  Some  improvement  might  have  been 
gained  had  additional  secondary  predictors  that  incorporate  past  information 
been  used  (e.g.,  the  past  12-hr  thickness  of  the  layer  between  the  constant- 
pressure  surface  being  "predicted"  and  the  next  lower  constant-pressure 
surface). 

Tests  of  the  Navy's  equations  on  independent  data  showed  that  the  best 
results  were  obtained  in  the  shorter  extrapolation  intervals  (100  to  50  mb 
and  50  to  30  mb),  and  we  use  these  equations  for  specifying  initial-guess  fields 
for  height  and  temperature  at  50  and  30  mb.  The  equations  did  not,  however, 
yield  results  as  satisfactory  in  the  longer  extrapolation  intervals  (200  to  100  mb, 
200  to  50  mb,  and  100  to  30  mb). 

The  hypothesis  to  be  tested  here  is  that  the  systematic  incorporation  of 
information  contained  in  present  and  recent-past  observations  at  mid-  and 
upper-tropospheric  levels  (as  well  as  that  contained  in  stratospheric  observa¬ 
tions)  will  yield  both  a  first-guess  field  and  an  analysis  at  stratospheric  levels 
which  are  significantly  superior  to  analyses  incorporating  only  information 
contained  in  stratospheric  observations  of  the  present  and  recent  past. 

The  analysis  technique  applied  to  the  initial-guess  fields  is  a  modifica¬ 
tion  of  the  technique  described  by  Cressman  [2] .  Limited  developmental  testing 
of  the  analysis  technique  was  performed  to  determine  near-optimum  values  for 
the  variables  used  in  the  analysis. 
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2.0  METHODS  FOR  COMPUTATION  AND  ANALYSIS 


2,1  Vertical-extrapolation  Screening-regression  Experiments 

A  program  that  derives  multiple  linear-regression  equations  through  the 
use  of  a  screening  procedure  has  been  written  by  Enger  and  Rodante  [3] .  The 
screening  procedure  selects,  from  a  set  of  possible  predictors,  a  subset  that 
contributes  significantly  and  independently  to  reducing  the  variance  of  the 
predictand.  The  predictor  having  the  highest  linear  correlation  with  the  specified 
predictand  is  selected  from  an  array  of  possible  predictors.  Next,  the  partial 
correlation  coefficients  between  each  of  the  remaining  predictors  and  the  pre¬ 
dictand  (holding  the  first  predictor  constant)  are  examined,  and  the  predictor 
associated  with  the  highest  coefficient  is  then  selected  as  the  second  predictor. 
Additional  predictors  are  selected  in  a  similar  manner.  This  procedure  is 
repeated  until  a  chosen  predictor  fails  to  explain  a  significant  additional  percent¬ 
age  of  the  remaining  variance  of  the  predictand  or  until  a  specified  maximum 
number  of  predictors  is  selected.  The  criterion  of  significance  employed  is  a 
modified  F-test  described  by  Miller  [5] . 

Table  2-1  is  a  list  of  the  predictands  and  their  respective  possible  predictors. 
In  the  table,  the  300-mb  height  and  temperature  each  appear  twice  as  predictands 
to  determine  whether  it  is  necessary  to  analyze  the  400-mb  level  to  obtain  a  good 
initial  guess  in  the  build-up  to  300  mb.  If  the  dependent-data  root-mean-square 
error  (RMSE)  in  extrapolating  from  500  to  300  mb  is  "sufficiently  less"  (arbitrarily 
determined)  than  the  combined  error  in  extrapolating  from  500  to  400  mb,  and 
from  400  to  300  mb,  then  one  could  extrapolate  from  500  to  300  mb  without 
sacrificing  the  accuracy  of  the  analysis  build-up. 

The  screening-regression  experiments  were  carried  out  (on  an  IBM  7090 
computer  at  United  Aircraft  Corporation  Research  Laboratories  in  East  Hartford) 
using  radiosonde  data  from  historical  Northern  Hemisphere  data  tapes  for  four 
midseason  months:  July  and  October  1958  and  January  and  April  1959.  The  data 
were  stratified  into  six  latitude  bands.  Appendix  A  lists  the  stations  used  in  each 
of  the  latitude  bands. 
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TABLE  2-1 

LIST  OF  "PREDICTANDS"  AND  POSSIBLE  PREDICTORS  USED  IN  VERT I  CAL -EXTRAPOLATION 

SCREEN  I NG- REGRESS  I  ON  EXPER I MENTS* 


Predictand 

Possible  predictors 

-pr 

O 

o 

NJ 

o 

4ooz_,2»  5oo-4ooh_|2>  500Z0,  500T0 

4ooTo 

400T_i2»  500-4oOT^^.12,  5OOZ0,  500Tq,  5oo-4ooHo 

300Zo 

300Z_,2,  500-300H_^2’  500Zq,  500Tq 

500Zq 

300Z_,2’  4oO-500H.^2»  ^OOZq,  400Tq 

500To 

300T_,2'  500-300T^^_12,  500Zq,  500Tq,  50O-5OOH0 

500To 

3OOT.12.  ^00-500T^^_^2»  *^00Zq,  400Tq,  4oO-300Ho 

25OZ0 

250Z_^2»  500-250H_,2'  ^OOZq,  300Tq 

25OT0 

250T_,2i  500-250T^^.12>  500Zq,  30OT0,  500-250Hq 

20OZ0 

200Z_^2»  25O-20OH.12.  25OZ0,  250Tq 

200Tq 

200T_,2»  250-200T^^_^2»  250Zq,  250Tq,  250-200Hq 

I5OZ0 

150Z.^2»  200-1  50H_^2»  200Zq,  200Tq 

I5OT0 

150T_^2»  200Zq,  200Tq,  200-1  50Hq 

IOOZq 

100Z_^2»  150-100H,^2»  1  50Zq,  1  50Tq 

IOOTq 

100T.^2»  ''50-100T^^.12,  1  50Zq,  150Tq,  150-100Hq 

*H,  T,  and  Z  are  thickness  (ft)  of  layer  between  constant-pressure 
surfaces,  temperature  (®C)  at  constant- pressure  surface,  and  height  (ft) 
of  constant- pressure  surface,  respectively.  Subscript  Q  designates 
observation  time,  subscript  -12  designates  12  hr  before  observation,  and 
subscript  m  designates  a  mean  over  a  layer. 
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The  regression  program  allows  for  several  predictands  and  their  predictors 
to  be  specified  at  one  time  by  generating  a  large  covariance  matrix  comprising 
all  predictands  and  predictors.  An  interchange  of  designation  of  variables  as 
predictors  or  predictands  for  different  e3q)eriments  is  also  possible  with  the 
generation  of  the  large  matrix.  The  ability  to  determine  several  regression 
equations  with  one  machine  run  results  in  a  marked  saving  of  computer  time  over 
deriving  each  equation  in  separate  runs. 

The  variability  of  times  at  which  a  station  reports  during  the  month  and 
the  variability  in  the  termination  level  of  the  sounding  (when  there  is  a  sounding) 
made  a  preprocessing  program  necessary  for  the  selection  of  cases  for  the 
predictands  [11] .  The  14  predictands  listed  in  Table  2-1  are  separated  into 
three  groups  of  six,  four,  and  four  (reading  from  top  to  bottom)  for  all  stations. 

If  a  station  is  to  be  selected  as  a  case  for  the  screening  program,  it  must  report 
all  parameters  in  its  group  (predictors  and  predictands);  if  even  one  is  missing, 
that  station  is  not  selected.  The  number  of  cases  for  each  group  and  each  latitude 
band  is  determined  for  the  month  in  question.  The  procedure  sometimes  excludes 
a  few  cases  that  would  be  included  were  individual  matrices  generated  for  each 
predictand,  but  the  advantage  in  saved  computer  time  far  outweighs  the  loss  of 
the  few  cases. 

At  the  other  extreme,  a  matrix  for  all  14  predictands  and  their  respective 
predictors  could  be  generated  -and  would  yield  a  further  saving  of  computer  time. 
This,  however,  would  require  that  all  14  predictands  and  their  predictors  be 
reported  and  would  result  in  the  loss  of  stations  whose  soundings  "topped  off" 
below  100  mb-as  many  as  118  cases  in  a  latitude  band.  Obviously,  this  would  be 
unsatisfactory.  The  breakdown  into  the  groups  of  six,  four,  and  four  provides  the 
most  feasible  solution  to  the  problem  of  selecting  enough  cases  while  minimizing 
computer  time. 
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2.2  Analysis  Procedure,  500-100  mb 


The  initial  guess  for  the  500-mb  height  analysis  is  the  previous  12-hr 
500-mb  height  analysis.  For  operational  purposes,  it  is  suggested  that  the 
initial-height  guess  at  500  mb  be  a  12-hr  dynamical  prognostic  chart  valid  for 
analysis  time.  The  initial  guess  for  the  temperature  field  at  500  mb  is  the 
previous  12-hr  temperature  analysis.  Once  again,  operationally,  if  a  12-hr 
dynamical  temperature  prognostic  chart  is  available,  it  should  be  used  as  the 
initial  guess. 

Station  observations  of  500-m.b  heights  and  winds  at  analysis  time  are 
used  successively  to  "correct"  the  initial-guess  field  heights,  using  a  technique 
very  similar  to  that  in  use  at  the  Numerical  Weather  Prediction  (NWP)  Unit  of 
the  National  Meteorological  Center  (NMC)  and  outlined  by  Cressman  [2] . 
Essentially,  the  technique  is  the  following. 

In  a  given  scan  through  the  data,  a  "correction"  is  computed  for  each  grid- 
point  within  a  specified  radius  R  of  any  reporting  station.  A  weighting  factor  is 
computed  for  each  station  as  a  function  of  the  distance  between  the  gridpoint 
being  corrected  and  the  station.  The  reported  wind  direction  and  speed  are  used 
in  conjunction  with  the  heights  to  help  establish  the  gradient  of  the  height  contours. 
[The  procedure  will  be  referred  to  here  as  the  successive-approximation  tech¬ 
nique  (SAT).]  The  "correction"  equation  [Eq.  (2-1),  below]  allows  for  the 
station  observation  to  be  weighted  relative  to  the  initial  guess  at  the  gridpoint  as 
well  as  to  the  winds,  and  hence  represents  a  change  in  the  equation  used  at  the 
NWP  Unit. 


W(d) 


Z  +  nZ  ^ 
x,y  obs 


mkfv. 


mkfu„ 


1  +  n  g 

Here,  W(d)  is  the  distance-weighting  function;  Z 


g 


(2-1) 


x.y 


is  the  initial-guess  value  of 


height  at  the  gridpoint  being  corrected;  is  the  observed  height  at  the  station; 
n  is  the  weighting  factor  of  the  station  observation  to  the  initial  guess  at  the  grid- 
point;  m  is  a  map-scale  factor;  k  represents  the  average  ratio  of  the  geostrophic 
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Fig.  2-1.  Schematic  representation  of  gridpoints  used  in 
curvilinear  interpolation. 


to  the  actual  wind  and  is  set  at  1.08  (the  value  suggested  in  Cressman's  paper 

[2] ;  g  is  the  acceleration  of  gravity;  u^  and  v^  are  the  components  of  the  wind 

along  the  i-  (horizontal)  and  j-  (vertical)  directions  of  the  grid,  respectively; 

x'  and  y'  are  the  distances  in  grid  units  in  the  i-  and  j-directions,  respectively, 

with  the  origin  taken  at  the  reporting  station;  and  S(x  ,y  )  is  obtained  by  fitting 

s  s 

a  curvilinear  surface  exactly  to  the  initial-guess  values  G(x,y)  at  the  four  grid- 
points  surrounding  the  observation  station. 

The  equation  for  a  curvilinear  surface  is  given  by 

S  =  a  +  bx  +  cy  +  dxy,  (2-2) 

where  x  and  y  are  distances  in  units  of  grid  intervals  in  the  i-  and  j-directions 
on  the  NWP  grid,  and 
a  =  G(0,0), 
b  =  G(1,0)  -  G(0,0), 
c  =  G(0,1).  -  G(0,0), 
d  =  G(l,l)  -  G(0,1)  -  G(1,0)  +  G(0,0), 

where  the  G's  are  as  indicated  in  Fig.  2-1.  The  total  "correction"  sqiplied  to  a 
gridpoint  is  simply  the  mean  of  the  corrections  computed  during  a  given  pass 
through  the  data. 

The  temperature  analysis  also  uses  SAT  in  arriving  at  a  final  field  of 
temperature  at  gridpoints.  For  temperature,  the  correction  equation  is 
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(2-3) 


C  =  W(d)  [T  -  S(x  y  )), 

a  9  a 

where  W(d)  and  S(x  ,y  )  are  as  in  Eq.  (2-1),  and  T  is  the  observed  value  at  the 
s  s  s 

station. 

The  500-mb  height  and  temperature  fields  thus  analyzed  form  the  base  for 
extrapolating  to  successively  higher  levels. 

For  analysis  of  height  and  temperature  at  the  next-higher  level  (400  or 
300  mb)  the  heights  are  specified  first,  using  the  results  of  the  screening- 
regression  experiments.  The  specified  400-  (or  300-)  mb  height  field  is  the 
initial  guess  for  the  400-  (or  300-)  mb  analyses  to  be  carried  out  in  the  same 
manner  as  that  described  above  for  the  500-mb  height  analysis.  It  appears  that 
analyzing  heights  before  temperatures  at  the  next-higher  level  is  more  useful 
than  analyzing  temperatures  first  because  the  winds  aid  in  arriving  at  a  height 
analysis,  and  the  resulting  height  analysis  may  then  be  used  to  compute  a 
current  thickness  field.  This  field,  in  turn,  is  very  likely  to  be  helpful  in  obtain¬ 
ing  a  more  accurate  temperature  analysis. 

From  400  (or  300)  mb,  the  analysis  proceeds  to  100  mb  in  steps  of  50  mb 
(but  bypassing  350),  with  the  height  analysis  being  performed  first  and  the  temper¬ 
ature  analysis  second.  The  results  of  the  regression  experiments  supply  the  initial 
guess  at  the  next-higher  level  for  the  heights  and  temperatures  at  gridpoints. 

Beginning  at  100  mb,  solar  radiation  errors  become  important.  The  100-mb 
height  and  temperature  data  used  in  the  analysis  are  corrected  for  radiation 
effects.  The  corrections  are  based  on  the  work  of  the  Stratosphere  Project  of  the 
U.S.  Weather  Bureau  under  the  direction  of  S.  Teweles  [12]. 

For  the  wind  analyses,  gradient  winds  are  computed  from  the  height  field 
at  every  level  and  are  used  as  an  initial  guess  to  a  direct  wind  analysis  using 
SAT.  Winds  are  resolved  into  u-  (east-west)  and  v-  (north-south)  components  for 
the  purposes  of  analysis,  but  the  capability  for  printing  wind  direction  and  wind 
speed  from  the  final  passes  of  the  u-  and  v-component  analyses  is  included  in  the 
analysis  program. 
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2.3  Analysis  Procedure,  100-30  mb 

From  100  mb,  the  Navy  vertical-extrapolation  regression  equation  is  used 
to  obtain  the  initial  guesses  of  height  and  temperature  at  50  mb.  These  50-mb 
fields  are  analyzed  in  the  same  manner  as  that  described  in  Section  2.2  (SAT)  for 
the  other  constant-pressure  surfaces.  The  50-mb  analyses  are  extrapolated  to 
30  mb  through  use  of  the  50-30-mb  Navy  extrapolation  coefficients  for  height  and 
temperature,  and  the  30-mb  analyses  use  SAT.  The  radiosonde  data  at  50  and  30 
mb  are  corrected  for  radiation  errors  as  discussed  above  for  the  100-mb  data. 
As  before,  gradient  winds  are  computed  from  the  height  analyses  and  are  used  as 
initial  guesses  to  direct  wind  analyses. 
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3.0  EXPERIMENTS  AND  EVALUATION  PROCEDURES 


3.1  Evaluation  of  Vertical-extrapolation  Regression  Equations,  with  Dependent 

and  Independent  Data 

The  regression  equations  developed  for  the  four  midseason  months  by  the 
program  discussed  in  Section  2.1  were  tested  on  independent  data.  (Radiation 
corrections  were  applied  to  the  dependent  and  independent  100-mb  data  where 
necessary.)  The  independent  data  were  selected  from  radiosonde  observations 
in  the  historical  Northern  Hemisphere  data  books  only  from  midseason  months  of 
years  not  included  in  the  developmental  data  sample.  Among  the  independent  data 
were  stations  not  included  in  the  original  dependent  sample.  This  allowed  the 
data  to  be  as  truly  independent  as  possible  with  respect  to  both  time  and  geography 
within  a  latitude  band.  Fifty  cases  were  chosen  as  independent  data  for  each  of 
the  six  latitude  bands  for  each  of  the  four  midseason  months. 

The  equations  were  tested  on  independent  data  for  the  months  for  which 
they  were  developed.  They  were  also  tested  on  data  from  the  months  on  each 
side  of  the  midseason  months  to  determine  whether  they  were  stable  and  useful 
for  those  months  as  well. 

3.2  Comparison  of  TRC  Equations  and  Navy  Equations 

As  suggested  in  Section  1.0,  it  has  been  our  belief  that  the  Navy  vertical- 
extrapolation  regression  equations  could  be  improved  with  the  use  of  additional 
predictors.  An  ejqjeriment  was  performed  to  test  this  hypothesis. 

In  the  experiment,  a  200-100-mb  height  regression-extrapolation  equation 
was  derived  for  the  25-35°  latitude  band.  This  extrapolation  interval  and  this 
latitude  band  were  selected  for  two  reasons;  (a)  the  200-100-mb  is  the  only 
extrapolation  interval  directly  comparable  to  the  intervals  for  which  we  are 
deriving  equations,  and  (b)  the  tropc^ause  is  most  likely  located  in  the  200-100-mb 
layer  in  January,  in  the  25-35°  latitude  band,  and  it  is  through  the  tropopause  region 
that  the  "prediction"  of  heights  is  most  difficult.  At  higher  latitudes,  the  tropopause 
is,  of  course,  usually  found  at  a  layer  lower  than  the  200-100-mb  layer. 


The  possible  "predictors"  specified  for  the  TRC  equation  are  200Zq, 

200Tp,  200-150H  150-100H  and  lOOZ  (The  symbols  and  subscripts 

are  identical  to  those  in  Table  2-1.)  The  Navy  regression  coefficients  were 
available  for  the  20-30*  and  30-40°  latitude  bands,  using  as  "predictors"  200Z^ 
and  200Tq. 

Both  equations  were  tested  on  the  same  independent  data  (100  cases)  for  the 
25-35*  latitude  band.  It  was  necessary  when  using  the  Navy  equation  to  use  two 
sets  of  coefficients,  one  for  the  stations  in  the  20-30*  latitude  band  and  the  other 
for  stations  in  the  30-40*  latitude  band.  The  TRC  equation  used  one  set  of 
coefficients  derived  for  25-35*  latitude  band. 

3.3  Comparison  of  Objective  and  Subjective  Wind  Analyses  at  100  mb 

Before  electronic  computers  became  available,  constant-pressure  charts 

were  necessarily  drawn  by  hand.  It  is  natural,  therefore,  to  compare  an  objective 
analysis  with  a  conventional,  subjective  hand  analysis  made  from  the  same  data. 

For  this  comparison,  objective  analyses  of  wind  at  100  mb  were  performed 
with  gradient  winds  computed  from  height  analyses  as  initial  guesses  for  three 
observation  times  in  January  1959.  The  direct  wind  analysis  was  carried  out  by 
SAT,  with  the  actual  wind  observations.  A  percentage  of  the  observations  was 
withheld  from  the  analysis  for  verification  purposes.  The  subjective  hand  analyses 
were  drawn  from  the  same  data  used  to  verify  both  kinds  of  analysis. 

In  determining  an  analysis  value  at  a  withheld  station,  a  curvilinear  surface 
[Eq.  (2-2)]  was  fitted  exactly  to  the  values  at  the  four  gridpoints  surrounding  the 
station.  It  was  necessary,  for  analysis  and  verification  purposes,  to  separate  the 
wind  into  its  two  components,  u  and  v.  (In  this  particular  case,  u  and  v  are  with 
respect  to  the  i-  and  j -directions,  respectively,  of  the  NWP  grid.) 

3.4  Developmental  Testing  of  the  Objective  Analyses,  and  Verification  by  the 

Areal-mean-error  Method 

In  the  objective  analysis  procedure,  several  parameters  may  be  varied  during 
any  given  analysis  (e.g.,  the  influence  radius  R  from  a  gridpoint,  and  the  weighting 
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factor  n,  which  weights  the  station  observation  of  height  relative  both  to  the 
initial  guess  at  a  gridpoint  and  to  the  winds).  Developmental  testing  consists 
in  assigning  several  values  to  these  parameters  and  determining  those  that  will 
minimize  the  analysis  errors.  The  analyses  are  verified  by  the  areal- mean- 
error  method  [14] .  This  method  s^proximates  an  area  integral  of  the  error 
over  the  analysis  area  by  withholding  some  percentage  of  the  observations  and 
verifying  back  to  the  "withheld"-  and  "analysis'-station  observations,  which 
are  equally  weighted. 

The  areal  root-mean-square  error  (ARMSE)  is  given  by 


fr-  2  -1  r- 

y  (4>  -  s)  p  y  (*  -  s 

ARMSE  =  - -  - - ; - —  +  ~ — - - - 

2l  p-^  27  p-^ 

L  w  w  ~a  a 


,2  -111/2 
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where  the  subscript  w  refers  to  quantities  at  the  location  of  observations  not 
used  in  the  analysis  computations  (i.e.,  the  withheld-station  observations),  4>  is 
an  observed  value,  s  is  an  analysis  value  obtained  by  fitting  a  curvilinear  surface 
exactly  to  the  values  at  the  four  gridpoints  surrounding  the  observing  station,  and 
p  is  the  density  of  analysis  stations  in  the  region  about  a  station.  Analysis  stations 
are  defined  as  stations  used  in  the  analysis  computations,  and  quantities  pertaining 
to  them  are  denoted  by  a  subscript  a. 

The  withheld-station  observations  approximate  the  maximum-error  points 
in  the  analysis  field,  but  the  number  of  withheld  stations  required  by  the  areal- 
mean-error  method  does  not  usually  equal  the  number  of  maximum-error  points 
in  the  analysis-error  field.  Therefore,  the  withheld  stations  must  be  selected  in 
a  manner  that  will  allow  them  to  represent  equally  the  entire  spectrum  of  maximum 
errors  over  the  analysis  area.  The  assumption  here  is  that  the  maximum  error  is 
inversely  proportional  to  the  density  of  the  data.  For  the  purpose  of  developmental 
tests,  in  which  a  series  of  maps  is  analyzed  in  one  computer  run,  it  is  advisable 
to  use  a  permanent  set  of  withheld  stations  to  obtain  an  unbiased  sampling  of  the 
maximum  errors  [14]. 


To  choose  the  permanent  set  of  withheld  stations,  a  measure  of  station 
density  was  computed  for  each  observation  for  a  series  of  20  observation  times. 

The  stations  were  then  ranked  in  order  of  increasing  density,  and  the  ranked  set 
of  stations  was  divided  into  groiqjs  of  K  stations  each.  One  station  was  selected 
from  each  groiq)  on  the  basis  of  the  frequency  with  which  it  reported  during  the 
series.  The  station  that  reported  most  frequently  throughout  the  series  was 
usually  selected,  provided  that  it  was  not  in  the  immediate  area  of  some  previously 
selected  station.  If  it  was,  then  the  station  that  reported  next  most  frequently  was 
chosen. 

Most  of  the  developmental  testing  was  done  on  a  North  American  grid. 

Limited  testing  was  done  on  a  hemispheric  grid. 

3.5  Comparison  of  "Build-up"  and  Persistence  as  Initial  Guess  to  100-mb  Analyses 

Our  major  aim  of  specifying  the  best  possible  initial  guesses  for  analyses  at 
100  mb  required  that  our  procedure  for  build-up  analysis  (see  Section  2.2)  be  com¬ 
pared  with  persistence.  Our  procedure  uses  vertical-extrapolation  regression 
equations  to  generate  initial  guesses  in  steps  to  100  mb.  An  e}q)eriment  was  con¬ 
ducted  to  evaluate  the  relative  accuracy  of  100-mb  temperature  and  wind  analyses 
using,  on  the  one  hand,  12-hr  persistence  as  a  first  guess  and,  on  the  other,  the 
build-up  first-guess  field.  In  this  experiment,  a  series  of  maps  was  analyzed  for 
five  consecutive  upper-air  observation  times.  The  initial  analysis  time  was 
1200Z  January  1,  1959.  Identical  sets  of  analysis  and  verification  observations 
were  used  for  the  analysis  for  each  type  of  first-guess  field.  The  analysis  area 
for  this  e3q)eriment  was  the  North  American  grid,  as  illustrated  in  Fig.  4-10. 

The  results  of  the  experiment  are  given  in  Section  4.0. 
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4.0  RESULTS  OF  EXPERIMENTS 


4.1  Evaluation  of  Vertical-extrapolation  Regression  Equations,  with  Dependent 

and  Independent  Data 

4.1.1  Midseason  Months 

After  each  predictor  had  been  selected  according  to  the  screening  procedure 
described  in  Section  2.1,  the  appropriate  regression  equation  was  derived.  This 
resulted  in  1512  equations,  (i.e.,  four  equations  for  each  height  in  each  latitude 
band  and  month,  and  five  equations  for  each  temperature  in  each  latitude  band  and 
month). 

There  remained  the  problem  of  selecting  the  most  useful  equation  for  each 
predictand.  The  criterion  for  selection  in  most  cases  was  the  modified  F-test 
for  significance  of  a  "predictor,"  described  by  Miller  [5] .  It  was  subjectively 
decided  to  eliminate  some  of  the  'predictors "-even  though  the  F-test  showed  them 
to  be  significant-whenever  the  developmental  sample  results  indicated  that  there 
would  be  little  to  gain  in  the  operational  use  of  an  additional  predictor. 

The  question  posed  in  Section  2.1  regarding  the  necessity  of  using  one  set 
of  vertical-extrapolation  equations  from  500  to  400  mb  and  another  set  from  400 
to  300  mb  rather  than  one  set  of  equations  from  500  to  300  mb  is  answered  by  the 
dependent  and  independent  data  results.  The  results  indicate  that  for  July,  the 
extrapolation  from  500  to  300  mb  may  be  performed  without  generating  unaccept¬ 
able  errors  in  the  analysis  build-up  (i.e.,  the  largest  300-mb  RMSEs  among  all 
latitude  bands  are  1.5rc  for  temperature  and  66.9  ft  for  height).  For  other  mid¬ 
season  months,  it  would  profit  the  analysis  to  include  the  400-mb  level  in  the  build¬ 
up  procedure.  If  it  is  deemed  necessary  to  extrapolate  to  400  mb  for  one  latitude 
band,  the  program  makes  it  mandatory  to  extrapolate  to  400  mb  for  all  latitude 
bands.  Hence,  there  are  312  vertical-extrapolation  regression  equations  for  the 
four  midseason  months  for  extrapolation  from  500  to  100  mb.  These  equations,  as 
well  as  the  Navy  equations  for  extrapolation  from  100  to  30  mb,  are  given  in 
Appendix  B. 
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The  results  on  dependent  and  independent  data  appear  in  graphical  form 
in  Figs.  4-1  through  4-8  and  are  presented  in  tabular  form  in  Appendix  C. 

In  general,  the  independent-data  tests  indicate  that  the  equations  are 
stable,  fur  the  independent  errors  are  close  to  the  values  of  the  dependent 
data  errors.  The  largest  errors  in  temperature  occur,  as  might  be  e:q)ected, 
in  the  layers  in  which  the  tropopause  is  frequently  found.  This  is  strikingly 
revealed  in  the  January  temperature  graphs  (Fig.  4-1),  which  show  dependent- 
data  temperature  errors  to  be  of  the  order  of  2.5°C  for  the  stations  between 
latitudes  35  and  55°  in  the  250-200-mb  layer.  Synoptic  e3q)erience  agrees  with 
these  results,  since  the  tropopause  is  often  found  in  the  250-200-mb  layer 
between  latitudes  35  and  55°  during  January. 

The  graphical  temperature  results  for  July  (Fig.  4-3)  show  that  the  largest 
errors  in  temperature  are  in  the  vicinity  of  the  150-100-mb  layer  for  latitudes 
0-45“N,  where  synoptic  ejqierience  places  the  tropopause  in  the  summer.  The 
error  curve  for  the  35-45°  latitude  band  suggests  that  it  is  also  often  found  in 
the  200-150-mb  layer  as  well.  In  the  more  tropical  portions  of  the  0-45°  band, 
the  tropopause  is  probably  frequently  above  100  mb,  and  the  results  indicated  in 
Lea's  paper  [4J  for  the  100-50-mb  temperature  equations  bear  this  out. 

Although  the  errors  are  largest  for  the  extr£q}olation  intervals  in  which 
the  trq}q>ause  is  frequently  located,  they  would  be  larger  if  past  12-hr  informa¬ 
tion  were  not  included  as  possible  "predictors."  This  statement  may  be  confirmed 
by  examining  the  results  of  the  comparative-accuracy  test  between  the  Navy  and 
TRC  equations  (see  Section  4.2).  This  test  demonstrates  that  the  inclusion  of 
additional  parameters  significantly  improves  the  earlier  Navy  equations.  Further 
confirmation  may  be  obtained  in  the  tabular  results  appearing  in  Appendix  C, 
which  include  a  listing  of  the  selected  predictors  (i.e.,  those  predictors  that  are 
statistically  significant  as  determined  by  the  modified  F-test  mentioned  in 
Section  4.1.1).  It  is  seen  that  for  the  "more  difficult"  extrapolation  intervals 
(those  in  which  the  tropopause  is  often  found),  the  past  thickness  and  past  height 
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Fig.  4-7.  Dependent  and  independent  data  results  with  TRC  vertical-extrapolation  regression  equations  for  height  for 


are  selected  for  the  height  extrapolation,  and  the  past  temperature  is  found  to 
be  significant  in  specifying  the  temperature. 

The  graphs  of  the  dependent-  and  independent-data  RMSE  for  height 
(Figs.  4-5  through  4-8)  indicate  an  {^parent  double  maximum  of  errors,  which 
is  most  noticeable  on  the  January  and  October  graphs. 

The  maximum  in  middle  latitudes  at  the  200-  and  250-mb  levels  may  be 
e}q)lained  by  the  frequent  presence  of  the  tropopause  in  those  regions.  The 
second  area  of  large  errors  is  at  the  100-mb  level  and  is  not  confined  to  one 
latitude  band  or  month.  Instead,  a  rather  flat  gradient  of  errors  is  evident 
across  all  latitude  bands  and  most  months  because,  despite  the  radiation  correc¬ 
tions  applied  to  the  data  at  100  mb,  instrumental  errors  persist.  These  errors 
are  of  the  order  of  100  to  200  ft  at  100  mb  [1] .  Temperature  observations  in  the 
stratosphere,  on  the  other  hand,  are  considered  to  be  accurate  to  within  1  or  2°C; 
hence  a  secondary  maximum  of  errors  across  all  latitude  bands  does  not  appear 
in  the  temperature  results. 

4.1.2  "Side"  Months 

The  equations  derived  for  the  midseason  months  were  tested  on  the  month 
on  each  side  of  each  midseason  month  to  determine  whether  the  equations  were 
stable  and  useful  for  those  months.  Twenty-five  cases  were  used  for  each  "side" 
month  for  testing  and  evaluation.  The  independent-data  RMSEs  for  the  predictands 
appear  in  tabular  form  in  Appendix  C  alongside  the  results  for  the  midseason 
months. 

The  evaluation  of  the  RMSEs  for  the  "side"  months  revealed  that  in  no 
case  does  it  appear  that  new  equations  need  be  derived  for  an  entire  latitude 
band.  Figure  4-9  depicts  the  "predictands"  in  each  latitude  band  and  month  for 
which  it  appears  that  it  is  necessary  to  derive  a  new  equation.  It  also  shows  17 
"predictands"  whose  RMSEs  indicated  that  there  was  some  question  as  to  whether 
a  new  equation  need  be  derived.  Additional  cases  are  necessary  to  determine  the 
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Fig.  4-9(b).  Predictands  for  which  new  equations  appear  to  be  necessary. 


percentage  of  these  "questionable  predictands"  that  will  need  new  equations. 

Some  further  testing  with  a  larger  sample  will  be  carried  out  for  a  portion  of 
the  'predictands"  that  appear  to  need  new  equations  as  well  as  for  all  of  the 
"questionable  predictands."  In  addition,  it  is  planned  to  test  on  a  larger  sample 
some  of  the  predictands  whose  RMSEs  indicate  that  new  equations  are 
necessary.  These  additional  tests  are  aimed  at  more  accurately  determining 
the  stability  of  the  equations  with  the  "side"  months. 

It  is  seen  that  there  are  24  'predictands"  for  which  new  equations  will 
probably  be  derived,  of  which  14  are  for  temperature  and  10  are  for  height. 

Of  the  10  probable  new  height  equations,  five  are  for  the  500-300-mb  extrap¬ 
olation  interval  (which  will  not  be  used  in  the  analysis  build-up  except  during 
June,  July,  and  August),  but  the  five  equations  are  e^qiected  to  be  used  in  a  prediction 
technique  that  is  based  on  the  regression  equations. 

Inasmuch  as  each  of  84  equations  for  each  midseason  month  (except  60  for 
July)  was  tested  on  two  months  other  than  the  midseason  month  for  which  it  was 
developed,  it  was  possible  that  624  new  equations  would  have  had  to  be  derived 
if  the  RMSEs  proved  too  great.  But  it  now  seems  necessary  to  derive  only  20  to 
40  new  equations  (about  5%  of  all  that  is  "possible"),  and  this  would  indicate  that 
the  equations  are  apparently  stable,  for  the  most  part,  for  a  season.  This  is 
especially  true  of  the  January  (winter)  equations  (Fig.  4-9),  which  were  used  on 
data  from  December  and  February.  For  December,  there  are  no  'predictands" 
that  appear  to  need  new  equations,  but  one  that  is  "questionable";  for  February, 
there  are  only  three  needing  new  equations,  and  three  that  are  "questionable." 

Two-thirds  (18  of  24)  of  the  equations  that  need  to  be  rederived  pertain  to 
the  spring  (17  equations)  and  fall  (18  equations),  which  are  recognized  climatol- 
ogically  as  transitional  seasons. 

4.2  Comparison  of  TRC  Equations  and  Navy  Equations 

Section  3.2  described  e^qieriments  for  comparing  the  Navy  regression 
equations  in  the  200-100-mb  extrapolation  interval  with  an  equation  derived 
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through  the  TRC  screening-regression  program.  Both  equations  were  tested  on 
the  same  independent  data. 

Of  the  five  possible  "predictors"  designated  at  TRC  (2002^,  200X^1 
200-150H  ,  150-100H  ,  and  1002  ),  all  except  150-100H  significantly 

reduced  the  variance  of  the  predictand  (1002^).  The  predictors  for  the  Navy 
equation  are  2002^  and  200Tp.  The  results  are  presented  in  Table  4-1. 

To  test  the  significance  of  the  improvement  shown  by  the  TRC  equation  in 
specifying  the  100-mb  height,  "Student's"  t-test  for  paired  comparisons  (described 
by  Wadsworth  and  Bryan  [15] )  was  performed  on  the  100  cases.  The  results  of 
this  test  indicated  that  the  difference  between  the  TRC  and  Navy  RMSEs  is 
significant  at  less  than  the  1%  level;  i.e.,  the  difference  is  highly  significant. 

In  the  100  cases,  the  TRC  "prediction"  was  better  than  the  Navy's  66  times,  and 
tied  with  the  Navy  "prediction"  eight  times. 

It  should  be  kept  in  mind  that,  in  our  regression  experiments,  we  specify 
shorter  extrapolation  intervals  than  does  the  Navy  (e.g.,  200-150  and  150-100  mb, 
rather  than  200-100  mb,  as  in  the  above  comparisons).  This  allows  an  analysis 
to  be  performed  at  the  150-mb  level  before  extrapolating  to  100  mb.  Use  of  the 
shorter  extrapolation  intervals  should  logically  result  in  lower  RMSEs  and,  in 
fact,  does  result  in  lower  RMSEs,  as  evidenced  by  the  results  for  1002^  in  the 
25—35°  latitude  band  for  January. 

4.3  Comparison  of  Objective  and  Subjective  Wind  Analyses  at  100  mb 

Table  4-2  presents  the  results  of  the  comparison  of  the  objective  and  sub¬ 
jective  wind  analyses.  It  is  evident  that4he  objective  analyses  could  generally 
be  considered  at  least  as  good  as  the  subjective  analyses.  The  analyses  themselves 
are  illustrated  in  Figs.  4-10  through  4-12.  The  observations  at  the  station  with¬ 
held  from  the  analyses  are  indicated  on  these  charts. 


29 


q: 

<  CD  t 

H  LU  -E 

<  q:  CD 

o  — 

I  Z  LU 

H  O  X 

f-  2  — 

I  LU  H  -O 
53-  O  <  E 
2  ^  I 
lU  UJ  O  Q 
^  CL  X  O 
CD  LU  «aC  «— 
<  Q  tr  I 
K  2  1“  Q 
-  X  O 
LU  CVJ 


5  10 


34 


Fig.  4-11.  Comparison  of  objective  and  subjective  analyses  of  vdnd  components  for 
OOOOZ,  13  Jan.  1959.  (a)  u  -component,  subjective  analysis. 


Fig.  4-11.  (c)  V  -component,  subjective  analysis. 


38 


Fig.  4-12.  Comparison  of  objective  and  subjective  analyses  of  wind  conqwnents  for 
1200Z,  13  Jan.  1959.  (a)  u  -conq)onent,  subjective  anal3rsis. 
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4.4  Developmental  Testing  of  the  Objective  Analyses,  and  Verification  by  the 

Areal-mean-error  Method 

With  major  concentration  placed  on  derivation  and  testing  of  the  vertical- 
extrapolation  regression  equations,  the  developmental  testing  of  the  analysis 
technique  was  necessarily  limited.  The  testing  nevertheless  has  enabled  us  to 
specify  for  the  variable  parameters  in  the  analysis  procedure  a  set  of  values 
that  will  produce  good- if  not  necessarily  optimum-analyses. 

The  three  products  at  each  constant-pressure  surface  treated  by  the 
build-up  analysis  procedure  (see  Sections  2.2  and  2.3)  are  height,  temperature, 
and  wind.  In  our  developmental  test  program,  the  major  emphasis  has  been 
placed  on  obtaining  the  best  wind  and  temperature  analyses,  whereas  the  height 
analyses  were  considered  primarily  as  a  means  for  generating  initial  guesses 
to  direct  wind  analyses;  i.e.,  the  assignment  of  different  values  to  the  variable 
parameters  in  the  height  analysis  has  been  directed  toward  minimizing  the  errors 
in  the  gradient  wind  field  as  computed  from  the  height  analysis  rather  than  toward 
minimizing  the  errors  in  the  height  field. 

The  following  sections  discuss  the  developmental  testing  and  the  verification 
statistics  obtained  from  the  areal-mean-error  method.  The  height  and  wind 
analyses  and  statistics  are  considered  separately  from  the  temperature  analysis 
and  statistics. 

4.4.1  Height  and  Wind  Analyses  and  Verification 

The  variable  par  imeters  in  a  height  analysis  are  the  influence  radius  R, 
a  weighting  factor  n,  the  "correction  type,"  and  smoothing. 

R,  the  infiuence  radius  from  a  gridpoint,  determines  the  stations  to  be  used 
in  £q}plying  a  correction  to  the  gridpoint  and  is  of  prime  importance  in  the  distance¬ 
weighting  factor 


W(d)  = 


R 


(4-1) 
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in  the  correction  equation.  W(d)  is,  therefore,  proportional  to  R.  From  develq)- 
mental  testing  performed  on  a  North  American  grid,  it  was  found  that  the  differences 
obtained  in  the  verification  statistics  by  varying  the  R's  within  reasonable  limits 
are  too  small  to  be  considered  of  any  importance.  Hence,  it  was  decided  to  use  the 
same  R's  used  by  the  NWP  Unit  in  their  500-mb  analyses-since  they  produce 
verifications  at  least  as  good  as  any  other  set  of  R's  tried.  The  radii  are  5.9,  3.6, 
2.2,  and  1.5.  Because  the  amount  of  available  data  diminishes  as  higher  levels 
are  analyzed,  the  NWP  analyses  above  500  mb  use  larger  radii  for  the  first  two 
phases  of  their  objective  hemispheric  analyses  (R  =  6.9  and  3.9  at  300  mb; 

ltd 

R  =  7.9  and  4.5  at  200  mb).  Although  it  may  be  necessary  to  use  larger  radii 

Ifd 

in  hemispheric  analyses  above  500  mb,  it  was  not  feasible  to  do  so  on  the  North 
American  grid  used  in  our  developmental  testing. 

n  is  a  weighting  factor  in  the  correction  equation  and  weights  the  station 
observation  relative  both  to  the  initial  guess  at  the  gridpoint  and  to  the  winds 
[see  Eq.  (2-1)] .  As  n  0,  the  winds  are  effectively  weighted  as  1. 

Three  possible  types  of  correction  may  be  applied  to  a  gridpoint  during  a 
height  analysis.  Type  1  corrections  use  only  the  reported  station  height  itself. 

T5T)e  2  corrections  use  both  the  reported  height  and  the  reported  wind.  Type  3 
corrections  use  only  the  wind. 

The  analysis  procedure  allows  for  a  choice  of  correction  on  each  pass. 
Although  it  would  appear  useful  to  include  all  available  data  for  all  scans,  experi¬ 
ence  has  shown  that  using  the  winds  during  the  first  pass  assists  in  obtaining 
a  specification  of  the  height  field  more  quickly  than  if  both  heights  and  winds  are 
used.  To  approach  the  correct  value  with  one  scan  through  the  data,  the  weighting 
factor  n  is  assigned  a  high  value  relative  to  the  initial-guess  value  at  all  levels 
in  the  analysis  build-up.  After  the  first  pass,  all  available  data  are  used,  and  the 
type  of  correction  applied  is  determined  solely  from  the  information  reported  by 
a  station. 
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After  the  fir^t  pass  of  the  analysis,  a  value  for  the  weighting  factor  n  must 
be  determined  which  will  result  in  a  near-q)timum  analysis  at  each  constant- 
pressure  surface.  (Since  all  available  data  are  used  for  passes  2  through  4,  n 
weights  the  reported  heights  relative  to  the  winds  as  well  as  to  the  initial-guess 
value  of  height  at  the  gridpoints.) 

It  has  been  concluded  [2,  13]  that  the  wind  information  is  a  better  determinant 
of  height  gradient  than  the  heights  themselves.  [The  winds  in  the  correction  equation 
(2-1)  for  height  analysis  are  geostrophic.]  The  NWP  correction  equation  originally 
[2]  weighted  the  winds  in  a  500-mb  analysis  4:1  over  the  heights.  This  has  since 
[9]  been  changed  to  8:1.  S.  Teweles  has  recommended  that,  in  objectively  analyzing 
100-mb  heights,  the  weighting  factor  of  winds  over  heights  should  be  increased  to 
32:1  and  suggests  [13]  that  the  factor  should,  perhaps,  be  64:1  at  50  mb.  He  reasons 
that  although  radiation  corrections  are  £q)plied  to  the  stratospheric  data,  instru¬ 
mental  errors  increase  as  the  radiosonde  balloon  rises.  The  winds,  therefore, 
should  be  a  more  reliable  indicator  of  the  height  gradient  than  the  heights  them¬ 
selves. 

Experiments  were  run  to  test  the  effect  of  varying  n  between  relatively  large 
limits  at  100  mb.  In  one  e3q>eriment  of  five  successive  observation  times  in 
January  1959,  the  weighting  factor  was  set  at  the  32:1  value  [n  =  0.031,  Eq.  (2-1)] 
suggested  by  Teweles.  In  another  experiment,  with  other  parameters  the  same, 
the  winds  were  weighted  only  four  times  more  heavily  than  the  heights  (n  =  0.25). 

The  results  are  presented  in  Table  4-3,  and  an  example  of  how  the  analyses  them¬ 
selves  compared  is  given  in  Fig.  4-13. 

InFig.  4- 13(b),  in  which  the  weighting  factor  of  winds  over  heights  is  4:1, 
the  trough  extending  southward  from  central  Canada  through  the  west  central  United 
States  is  deeper  than  the  trough  in  Fig.  4- 13(a),  in  which  the  weighting  factor  is 
32:1.  The  difference  in  height  of  the  trough  between  the  two  maps  is  more  than 
200  ft  in  the  northern  part,  more  than  100  ft  over  the  west  central  United  States, 
and  negligible  over  the  extreme  southwestern  United  States. 
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TABLE 

WINO-ANALYSIS  VERIFICATION  STATISTICS^ 
(a)  Root-naan-aqaura  arror,  u-co^ponant 


W6E 

RMSE-A 

RM5E-M  { 

Map 

n  «  0.23 

n  -  0.031 

n  -  0.23 

n  «  0.031 

n  -  0.25 

n  -  0.031 

1 

9.710 

10.965 

8.062 

8.337 

11.100 

13.105 

2 

11.523 

12.311 

9.890 

10.918 

12.9*19 

13.922 

3 

10.237 

10.654 

7.945 

8.938 

12.101 

12.126 

4 

11.319 

13.169 

9.429 

10.502 

12.932 

15.377 

5 

10.745 

13.019 

8.838 

10.112 

12.358 

15.383 

Ovar-al 1 

10.743 

12.113 

6.270 

6.9*19 

6.724 

9.922 

(b)  Root*me&n>squar0  arror»  v-cofliponant 


Map 

RM5E 

RMSE 

-A 

RHSE-W 

n  -  0.25 

n  »  0.031 

n  -  0.25 

n  -  0.031 

n  -  0.23 

n  *  0.031 

mm 

8.466 

8.481 

5.603 

5.872 

10.609 

10.457 

9.37*1 

9.700 

5.991 

6.347 

11.824 

12.159 

H 

13.937 

1*1.707 

6.692 

7.537 

18.536 

19.382 

HI 

9.337 

9.495 

8.270 

9.031 

10.291 

9.935 

mn 

9.361 

11.031 

8.6o4 

9.293 

10.058 

12.328 

10.337 

10.950 

5.110 

5.55*1 

8.986 

9.*157 

(c)  Averaga-vactor  wind  arror 


AVWE 

AVWE-A 

AVME-H 

n  -  0.25 

n  «  0.031 

n  -  0.23 

n  -  0.031 

n.-  0.25 

n  -  0.031 

12.896 

13.878 

6.957 

7.212 

10.859 

11.857 

14.855 

15.831 

8.177 

8.930 

12.401 

13.072 

B 

17.293 

18.160 

7.346 

8.268 

15.655 

16.169 

14.673 

16.236 

8.870 

9.796 

11.688 

12.947 

5 

14.250 

17.064 

8.722 

9.712 

11.269 

14.030 

Ovar*a| 1 

14.909 

16.329 

6.088 

8.896 

12.524 

13.693 

(d)  Avaraga  wind  spaad  on  final  pass 


Map 

AHS 

1 

38.527 

2 

40.271 

3 

42.535 

4 

43.391 

5 

45.074 

^Comparison  of  arrors  af  100  mb  for  winds  walghtad  ll:1  ovar  heights  (n  «  0.2^) 
against  arrors  for  winds  weighted  32:1  over  heights  (n  ■  0.031)  in  analysis  procedure. 
tAII  RMSESi  AVWEs,  and  AWSs  are  in  knots. 
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procedure,  (a)  Winds  weighted  32:1  over  heights. 


530  526  522  518  514  510506 


Note  also  that  the  ridge  over  eastern  and  northeastern  Canada  that  appears 
in  Fig.  4-13(b)  does  not  t^pear  in  Fig.  4-13(a).  This  ridge  is  ^parently  an  exten¬ 
sion  of  the  ridge  off  the  East  Coast.  The  minor  ridge  in  Fig.  4-13(a)  over  the 
western  edge  of  Hudson  Bay  is  not  evident  in  Fig.  4- 13(b).  The  reality  of  this 
ridge  may  be  questioned  on  the  basis  of  its  limited  areal  extent. 

The  results  of  the  combined  experiment,  as  evidenced  by  the  verification 
statistics  (Table  4-3)  using  the  areal-mean-error  method,  are  somewhat  surpris¬ 
ing,  especially  since  the  parameter  being  verified  is  the  gradient  wind  as  computed 
from  the  height  analysis.  In  general,  the  100-mb  analysis,  with  the  winds  weighted 
4:1  over  the  heights,  results  not  only  in  better  verification  statistics  but  also  in 
what  is  probably  a  more  realistic  analysis  than  the  analysis  with  the  winds  weighted 
32:1  over  the  heights. 

The  reported  winds  were  plotted  in  Fig.  4- 13(b)  to  ascertain  whether  cross¬ 
contour  flow  existed.  Teweles  [13]  found  it  in  100-mb  hemispheric  analyses  in 
which  the  winds  were  weighted  at  a  ratio  less  than  32:1  over  the  heights.  As  can 
be  seen,  there  is  very  little  cross-contour  flow  on  this  particular  North  American 
analysis. 

Apparently,  analyses  that  weight  the  winds  lightly  give  consistently  better 
results  than  analyses  that  weight  the  winds  heavily  because  the  radiation-corrected 
100-mb  heights,  though  containing  instrumental  errors,  are  able  to  specify  the 
contour-gradient  field  well  enough  to  yield  reasonably  good  computation  of  gradient 
winds.  It  is  well  known  [9]  that  the  winds  are  subgeostrophic  in  troughs  and  super- 
geostrophic  in  ridges;  hence,  the  gradient  winds  are  a  better  measure  of  the  true 
winds  in  troughs  and  ridges  than  are  the  geostrophic  winds. 

The  lighter  weighting  of  the  winds,  though  seemingly  a  more  desirable  feature 
in  a  100-mb  height  analysis  on  a  North  American  grid,  may  not  be  as  desirable  for 
analyses  of  hemispheric  charts.  The  data  are  relatively  dense  on  a  North  American 
grid,  and  the  radiosonde  instruments  are  more  compatible  over  North  America  than 
over  any  other  geographical  area.  Foreign  radiosondes  are  not  compatible,  and. 
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although  corrections  for  radiation  effects  have  been  develc^ed  for  them  by 
Rothenberg  and  Teweles  [8] ,  changes  in  instrumentation  (which  have  not,  to  date, 
been  accounted  for)  lead  to  a  lack  of  confidence  in  the  height  and  temperature  data 
reported  by  other  countries.  Therefore,  the  tendency  is  to  weight  the  winds  in  a 
hemispheric  100-mb  analysis  more  heavily  than  the  4:1  ratio  suggested  by  the 
verification  statistics  for  use  on  a  North  American  gprid. 

Table  4-4  gives  suggested  values  for  n  for  use  in  hemispheric  analyses  at 
various  constant-pressure  surfaces. 

Smoothing,  the  remaining  variable  parameter,  eliminates  the  small-scale 
"wiggles"  that  generally  appear  in  an  analysis.  Its  evaluation  here  is  particularly 
important  because  if  each  constant-pressure  surface  is  smoothed  after  the  final 
pass  through  the  data,  and  the  regression  equations  are  used  to  extrapolate  upward 
to  obtain  an  initial  guess  at  the  next-higher  level,  the  effect  of  smoothing  becomes 
cumulative  [7] .  How  much  effect  the  accumulated  smoothing  has  on  the  analyses 
at  the  levels  in  which  we  are  most  interested  (100,  50,  and  30  mb)  remains  to  be 
determined.  A  decision  as  to  what  type  of  smoothing  operator  to  employ,  and  when 
and  where  to  use  it,  must  be  made  on  the  basis  of  the  verification  statistics  obtained 
during  our  developmental  test  runs. 

The  NWP  Unit  vises  a  complex  9-point  operator  designed  by  Schuman  [10]. 
Their  experience  has  been  that  a  5-point  operator  has  imdesirable  characteristics. 
The  9-point  operator  effectively  eliminates  the  subsynoptic-scale  analysis  "wiggles" 
of  less  than  two  NWP  grid  intervals.  Further,  it  has  some  effect  in  reducing  the 
amplitude  of  waves  with  wavelength  between  two  and  four  grid  intervals  and  leaves 
waves  of  wavelength  greater  than  four  grid  intervals  virtually  vmchanged. 

Time  was  not  available  in  our  limited  developmental  testing  either  to  design 
a  smoothing  operator  such  as  the  NWP  Unit's  complex  operator  or  to  ad{q)t  their 
smoothing  program  to  our  analysis  program.  Hence,  it  was  necessary  to  design  a 
simplified  9-point  operator  of  the  form 

<Z)  =  +  k^(Z^  +  Z^  +  Zg  +  Zg)  +  +  Zg  +  Zg  +  Z^). 

(4-2) 
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TABLE  4-4 

WEIGHTING  FACTOR  n  USED  IN  BUILD-UP  ANALYSIS  PROCEDURE 


Leve  1 , 
mb 

Pass  1 , 

R  =  5.9 

Pass  3, 

R  =  2.2 

Pass  4, 

R  =  1 .5 

500 

5.0 

0.25 

0.25 

0.25 

4oo 

5.0 

0.25 

0.25 

0.25 

500 

5.0 

0.25 

0.25 

0.25 

250 

5.0 

0.25 

0.25 

0.25 

200 

5.0 

0.125 

0.125 

0.125 

150 

5.0 

0.125 

0.125 

0.125 

100 

5.0 

0.063 

0.063 

0.063 

50 

5.0 

0.031 

0.031 

0.031 

30 

5.0 

o.oei 

0.021 

0.021 

li 

I! 

II 
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Fig.  4-14.  The  9-point  mesh  used  for  smoothing, 
where  the  subscripts  refer  to  gridpoints  in  Fig.  4-14,  and  k  ,  k  ,  and  k  are 

X  A  V 

constants,  the  optimum  values  of  which  may  be  determined  through  testing  of 
various  combinations. 

The  main  difference  between  TRC 's  9-point  operator  and  that  used  by  the 
NWP  Unit  is  that  the  TRC  operator  uses  one  set  of  real  constants  and  one  scan 
over  the  grid.  This  results  in  the  central  point's  being  influenced  only  by  the 
eight  surrounding  points.  The  operator  designed  by  Schuman,on  the  other  hand, 
uses  three  elements  (indices),  two  of  which  are  conjugate  complex  indices,  and 
three  scans.  This  effectively  smooths  the  field  three  times,  and  all  of  48  surroxmd- 
ing  points  influence  a  smoothed  point. 

Four  e:q)eriments  using  the  build-up  analysis  procedure  were  run  to  test 
the  effects  of  smoothing  and  the  sensitivity  of  height  analyses  to  different  smooth¬ 
ing  constants  in  Eq.  (4,-2).  In  all  four  experiments,  the  other  variable  parameters  in 
the  analysis  procedure  were  kept  constant;  only  the  specifications  for  smoothing 
were  varied.  Temperature  and  wind  analyses  were  not  smoothed.  Each  of  the 
e^eriments  was  run  on  data  from  the  first  five  observation  times  of  January  1959, 
and  the  analyses  were  performed  on  a  North  American  grid.  In  all  cases  where 
smoothing  was  applied,  it  was  applied  to  the  last  pass  of  the  analysis. 

Experiment  1.  Height  analyses  at  all  constant-pressure  surfaces  were 
smoothed,  with  the  smoothing  constants  k  ,  k  ,  and  k  [Eq.  (4-2)]  set  equal  to 
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0.8>  0.035,  and  0.015,  respectively.  (For  smoothing  along  the  boundary  (but  not 
at  corner  points),  Fq.  (4-2)  was  suitably  modified  to  include  the  five  surrounding 
points,  and  smoothing  at  corner  points  was  accomplished  using  only  three  sur¬ 
rounding  points.] 

Experiment  2.  Smoothing,  using  the  same  values  of  k  ,  k  ,  and  k  specified 

X  A  U 

for  experiment  1,  was  carried  out  for  height  analyses  at  only  the  300-,  200-,  and 
100-mb  levels. 

Experiment  3.  No  smoothing  was  performed. 

Experiment  4.  Height  analyses  at  all  constant-pressure  surfaces  were 
smoothed,  with  k^  =  0.6,  k^  =  0.075,  and  k^  =  0.025. 

The  resulting  analyses  and  verification  statistics  indicated  the  following. 

(a)  Smoothing  at  all  levels  with  k  =  0.6,  k  =  0.075,  and  k  =  0.025  resulted 

X  ^  3 

in  too  heavy  a  smoothing,  filling  cyclones  and  troughs,  and  weakening  anticyclones 
and  ridges  to  an  undesirable  extent. 

(b)  No  smoothing  at  all  was  not  satisfactory,  as  evidenced  by  the  verification 
statistics  presented  in  Table  4-5  for  experiments  1,  2,  and  3.  The  initial  guesses 
(IG'  in  Table  4-5  refer  to  the  gradient  winds  computed  from  the  last  pass  of  the 
height  analysis,  and  finai  pass  (FP)  refers  to  the  last  pass  of  the  wind  analysis  that 
uses  wind  observations  to  correct  the  initial  guesses.  It  is  seen  that  the  over-all 
average-vector  wind  errors  are  consistently  larger  (above  500  mb)  for  the  no¬ 
smoothing  experiment  (3)  than  for  either  of  the  other  two  experiments.  Further 
evidence  supporting  the  conclusion  that  no  smoothing  is  unsatisfactory  is  given  by 
Figure  4-15,  which  compares  wind  speeds  resulting  from  a  height  analysis  that 
was  smoothed  after  the  final  pass  and  wind  speeds  resulting  from  a  height  analysis 
that  was  not  smoothed.  The  unsmoothed  analysis  appears  to  produce  unrealistic 
gradients,  resulting  in  wind  speeds  that,  at  isolated  gridpoints,  are  excessive  in 
comparison  with  those  at  nearby  gridpoints. 
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TABLE  4-5 

over-all  WINO-VERIFICATION  statistics  for  5  hr*t 


12.84  3-78  17-76  9.7°  2.4o  |  13.50  16.09  60.71 


from  unsmoothed  height  analsrsis. 
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(c)  The  verification  statistics  listed  in  Table  4-5  show  that  there  are 
negligible  differences  in  the  average-vector  wind  error  at  100,  50,  and  30  mb  for 
experiments  1  (smoothing  at  all  levels)  and  2  (smoothing  at  300,  200,  and  100  mb). 
These  statistics  fail  to  answer  the  question  of  whether  smoothing  at  all  levels  in 
a  build-up  analysis  is  preferable  to  smoothing  at  selected  levels.  The  analyses 
themselves  do  not  really  answer  this  question  either,  for  they  appear  to  be  quite 
similar  over  most  regions  of  the  analysis  area.  There  is  a  tendency,  however, 
for  a  100-mb  height  analysis  that  has  resulted  from  smoothing  at  selected  levels 
to  define  troughs  and  ridges  more  clearly  than  an  analysis  that  is  smoothed  at 
all  levels.  This  tendency  is  demonstrated  in  Fig.  4-16.  It  is,  therefore,  left  to 
the  eventual  user  to  determine  which  of  the  two  alternatives  he  prefers  for  his 
analysis  procedure.  There  does  not  seem  to  be  much  question  that  smoothing  is 
necessary  at  the  100-mb  level  and  is  also  desirable  at  least  at  other  levels  below 
100  mb.  Above  100  mb,  a  light  smoothing  should  be  used  at  30  mb,  since  the 
verification  statistics  in  Table  4-5  indicate  that  smoothing  results  in  somewhat 
smaller  errors. 


In  computing  gradient  winds  from  height  analyses,  a  problem  became  evident 
after  the  initial  computer  riui.  The  gradient  winds  computed  below  latitude  30°N 
resulted  in  excessively  high  west-wind  components  that  turned  out  to  be  unrealistic 
in  view  of  the  contour  spacing.  The  gradient-wind  equation  is 


V  = 

gr 


fr 

2 


4V 

1 

fr 


1/2 


1 


(4-3) 


where  f  is  the  Coriolis  parameter,  r  is  the  radius  of  curvature  of  the  contours  in 

feet  per  grid  unit,  and  V  is  the  geostrophic  wind  in  feet  per  second.  (The  Coriolis 

®  -5  -1 

parameter,  f,  is  ecjual  to  2u  sin  tp  =  14.584  x  lO  sin  (p  rad  sec  ,  where  <p  is  in 

degrees  of  latitude.)  Apparently,  the  excessively  high  west-wind  components  occur 

when  the  square-root  term  in  Eq.  (4-3)  has  a  negative  value,  in  which  case  the 

wind  is  confuted  from  -fr/2. 
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smoothing  was  performed  at  any  level  in  the  analysis  build-iq>. 


A  change  in  the  program  to  compute  geostrophic  winds  when  the  square- 
root  term  in  Eq.  (4-3)  is  negative  helped  somewhat,  but  fictitiously  high  west- 
wind  components  continued  to  appear  in  low  latitudes.  An  investigation  traced 
their  cause  to  the  occurrence  of  a  slightly  positive  square-root  term  in  Eq.  (4-3), 
resulting  in  computation  of  gradient  winds  primarily  from  the  -fr/2  term.  To 
solve  this  problem,  geostrophic  winds  were  computed  at  those  gri(^oints  where 
the  curvature  was  anticyclonic  (r  negative)  south  of  latitude  30‘’N.  Results  then 
indicated  no  fictitiously  high  west-wind  components,  nor  were  there  any  discon¬ 
tinuities  in  the  wind  field  across  the  30th  parallel.  The  procedure  of  computing 
geostrc^hic  winds  as  initial  guesses  when  the  contour  curvature  is  anticyclonic 
south  of  30'’N  was  incorporated  into  the  wind-analysis  program  and,  as  a  result, 
no  further  problems  were  encoimtered  in  the  developmental  testing  of  the  wind 
analysis. 

The  only  variable  parameter  in  our  direct  wind  analysis  is  R,  the  influence 
radius.  Initially,  four  passes  through  the  data  were  tried,  but  verification  statistics 
indicated  that  little,  if  any,  improvement  was  gained  from  the  fourth  pass.  Only 
limited  developmental  testing  of  the  direct  wind  analysis  was  performed.  Varying 
the  R's  for  the  u-  and  v-components  (within  reasonable  limits)  did  not  result  in 
any  marked  difference  in  the  verification  statistics,  hence  the  choice  of  a  series 
of  R's  is  rather  arbitrary.  However,  in  choosing  the  R's,  it  should  be  realized 
that  the  u-component  of  the  wind  generally  has  a  longer  wavelength  than  the 


TABLE  4-6 

R,  INFLUENCE  RADII,  USED  IN  DIRECT  WIND  ANALYSES 


Fig.  4-17.  100- mb  temperature  analyses,  (a)  No  smoothing  performed. 


v-component,  so  the  R's  for  the  u-component  should  be  somewhat  larger  than 
those  for  the  v-component. 

Table  4-6  contains  a  set  of  R's  that  has  given  satisfactory  results. 

4.4.2  Temperature  Analysis  and  Verification 

The  variable  parameters  in  the  temperature  analysis  are  R  (the  influence 
radius)  and  smoothing.  Since  the  station  density  is  exactly  the  same  for  the 
temperature  observations  as  it  is  for  the  height  observations,  the  same  R's  used 
by  the  NWP  Unit  for  their  500-mb  analyses,  being  as  satisfactory  as  any  other 
set  that  was  tried,  were  used  for  our  temperature  analysis. 

An  e>q)eriment  was  run  to  determine  the  value  of  smoothing  a  temperature 
analysis  at  100  mb.  The  verification  statistics  for  five  observation  times  showed 
little  difference  between  smoothed  and  unsmoothed  temperature  analyses:  the 
smoothed-analysis  RMSE  was  1.32°C,  and  the  unsmoothed-analysis  RMSE  was 
1.40°C.  The  actual  analyses  exhibited  negligible  differences,  as  may  be  seen 
from  Fig.  4-17.  It  is  recommended  that  the  temperature  field  not  be  smoothed 
for  a  build-up  analysis  procedure. 

4.5  Comparison  of  Build-up  and  Persistence  as  Initial  Guess  to  100-mb  Analyses 

Five  consecutive  sets  of  observations  were  analyzed  by  two  methods.  Each 
method  was  identical  in  every  respect,  with  the  exception  of  the  specification  of 
the  initial-guess  field.  That  is,  the  same  set  of  analysis  and  verification  stations 
and  the  same  analysis  parameters  (i.e.,  R,  number  of  passes,  n,  etc.)  were 
employed  in  both  instances.  The  initial-guess  field  for  the  first  set  of  analyses 
was  defined  as  the  analysis  of  observations  taken  12  hr  earlier.  For  the  second 
set  of  analyses,  the  first-guess  field  was  defined  as  that  derived  by  the  Explication 
of  the  appropriate  vertical-extrapolation  equations,  as  outlined  in  Section  2.0. 

The  verification  results  are  summarized  in  Table  4-7. 


It  can  be  seen  that,  for  both  wind  and  temperature,  the  measure  of  error  is 
smaller  at  each  stage  of  the  analysis  for  the  build-up  method.  To  test  the  significance 


TABLE  4-7 

VERIFICATION  STATISTICS  AT  100  mb 
FROM  PERSISTENCE  AND  BUILD-UP  ANALYSIS  PROCEDURES 


Over-all  vector-wind  error,  knots 

Temperature 

RMSE,  'C 

Pass 

Persistence 

Bui  Id- up 

Persistence 

Bui  Id-up 

0 

16.21 

12.70 

2.5 

1.9 

1 

15.25 

11 .00 

2.0 

1.7 

2 

11 .81 

10.15 

1.7 

1.5 

5 

11.54 

9-79 

1.5 

1.5 

4 

— 

— 

1.4 

1.2 

of  the  indicated  improvement,  "Student's"  t-test  for  paired  comparisons  [15J  was 
performed  on  the  five  cases.  The  results  of  this  test  indicated  that  the  difference 
in  error  between  the  two  analysis  methods  was  significant  at  less  than  the  5% 
level  for  the  wind  analysis  and  not  significant  at  the  5%  level  for  the  temperature 
analysis. 

This  e3q)eriment  should  indicate  the  results  to  be  expected  in  regions  of 
relatively  dense  data.  To  complete  the  evaluation,  it  is  necessary  to  conduct  a 
similar  experiment  in  which  sparse-data  regions  of  the  hemisphere  are  simulated. 
Such  an  e3q)eriment  has  been  attempted;  however,  the  results  were  inconclusive 
because  of  a  bias  introduced  by  errors  in  the  data  furnished  TRC. 


5.0  SUMMARY  AND  CONCLUSIONS 


Testing  of  vertical-extrapolation  regression  equations  derived  from  mid¬ 
season  radiosonde  data  on  independent  data  from  other  midseason  months  yielded 
excellent  results.  These  equations  should  be  useful  for  specifying  good  initial- 
guess  fields  of  height  and  temperature  for  a  build-up  objective-analysis  procedure. 

Tests  of  the  equations  with  data  from  the  months  on  each  side  of  each  mid¬ 
season  month  revealed  that  only  about  5%  of  the  equations  produced  errors  large 
enough  to  warrant  new  equations.  Two-thirds  of  these  5%  pertained  to  months  in 
spring  or  fall  (transitional  seasons).  The  vertical-extrapolation  equations  developed 
for  a  midseason  month  appear  to  be  characteristically  stable  for  the  other  months 
of  the  season. 

For  the  200-100-mb  extrapolation  interval  a  vertical-extrapolation  regres¬ 
sion  equation  that  uses  present  and  past  data  was  tested  against  an  equation 
(ieveloped  by  the  Navyjthat  uses  only  present  data.  The  equation  that  uses  present 
and  past  (12-hr)  data  gave  significantly  better  results  (as  determined  by  a 
"Student's"  t-test)  on  independent  data,  demonstrating  the  value  of  past  12-hr 
information. 

Limited  developmental  testing  of  the  build-up  objective-analysis  procedure 
(incorporating  the  vertical-extrapolation  equations)  indicated  the  following. 

(a)  Varying  the  set  of  R's  (the  influence  radii)  within  reasonable  limits 
in  any  of  the  three  objective  analyses  (height,  temperature,  wind)  did  not  result 
in  significantly  changed  verification  statistics.  The  set  of  R's  used  by  the  NWP 
Unit  in  their  500-mb  analyses  gave  results  as  good  as  any  other  set  tried. 

(b)  Wind  observations  are,  in  general,  more  reliable  than  height  observations, 
and  wind  observations  should  be  weighted  more  heavily  than  height  observations 

in  the  correction  equation  for  heights  at  gridpoints.  A  wind-to-height  weighting 
ratio  of  4:1  is  recommended  for  height  analyses  in  the  middle  troposphere,  and 
this  ratio  should  be  increased  to  at  least  16:1  at  100  mb.  The  ratio  at  100  mb  is 
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dependent  on  the  geographic  area  of  the  analysis  and  the  compatibility  of  radio¬ 
sonde  stations  in  that  area.  For  North  America,  where  the  data  are  dense  and 
the  compatibility  of  radiosonde  instruments  is  high,  a  4:1  ratio  was  found  to 
give  better  verification  statistics  than  a  32:1  ratio.  For  hemispheric  analyses, 
on  the  other  hand,  Teweles  [13] ,  after  considerable  experimentation,  suggested 
the  32'.1  ratio.  Incompatibility  of  radiosonde  instruments  and  the  variation  in  data 
density  over  the  entire  hemisphere  argue  for  the  use  of  a  32:1  ratio. 

(c)  If  the  objective  is  a  stratospheric  analysis  at  the  top  of  the  build-vg), 
heights  need  not  be  smoothed  at  all  constant-pressure  surfaces  but  only  at  the 
stratospheric  levels  and  at  two  (or  more)  levels  below  100  mb  (and  at  any  other 
level  that  may  be  of  particular  interest). 

(d)  Smoothing  does  not  appear  to  be  necessary  to  any  temperature  analysis. 

A  comparison  test  was  performed  on  a  North  American  grid  to  evaluate 
the  relative  accuracy  of  100-mb  wind  and  temperature  analyses  as  obtained  from 
two  different  types  of  initial-guess  fields.  One  initial-guess  field  was  given  by 
12-hr  persistence,  and  the  other  was  generated  by  the  build-up  analysis,  using 
vertical-extr£g)olation  equations.  Results  of  the  test,  as  determined  from  the 
areal-mean-error  verification  statistics,  indicated  that  the  build-up  procedure 
produced  wind  RMSEs  that  were  significantly  lower  than  those  given  by  using 
persistence  as  an  initial  guess.  "Student's"  t-test  for  paired  comparisons  was 
applied  to  the  differences  in  the  RMSEs  and  indicated  the  build-vg)  errors  to  be 
significantly  lower  at  the  5%  level.  For  temperature,  the  build-up  procedure 
produced  lower  errors  than  persistence,  but  the  difference  was  not  significant 
at  the  5%  level. 
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APPENDIX  A.  STATIONS  USED  IN  SCREENING-REGRESSION  EXPERIMENTS 


Latitude  0-25"  N 


Block 

Station 

Deg.  lat. 

Oeg.  long. 

72 

201 

24.55N 

3l .80W 

76 

644 

20.97N 

89.65W 

78 

118 

21 .47N 

71-15W 

78 

325 

23.15N 

82.35W 

78 

367 

19.9ON 

75.1 5W 

78 

397 

1  7.93N 

76.78W 

78 

501 

17.4ON 

85.93W 

78 

526 

18.45N 

66.  low 

78 

866 

18.03N 

63.  low 

78 

967 

10.68N 

61.62W 

78 

988 

12.I8N 

68.98W 

91 

115 

24.78N 

i4i.33E 

91 

217 

13.55N 

144.83E 

91 

245 

I9.28N 

166.65E 

91 

250 

11.55N 

162.33E 

91 

275 

I6.75N 

I69.52W 

91 

285 

I9.73N 

155.07W 

91 

4o8 

7.35N 

134.48E 

91 

413 

9.52N 

138.1 3E 

Latitude  25-35'’N 


Block 

Station 

Deg.  iat. 

Deg.  long. 

47 

931 

26.55N 

127. 75E 

60 

119 

34.5ON 

6.6OW 

62 

oil 

32.9ON 

13.28E 

72 

208 

52.9ON 

8O.OOW 

72 

226 

52.3ON 

86.40W 

72 

252 

28.98N 

89.57W 

72 

235 

32.33N 

90.22W 

72 

24o 

50.22N 

93.1 5W 

72 

250 

25.92N 

97.47W 

72 

259 

32.83N 

97.05W 

72 

261 

29.33N 

100.88W 

72 

270 

31  .80N 

106.40W 

72 

290 

52.73N 

117. 17W 

72 

54o 

54.73N 

92.23W 

74 

794 

28.47N 

80. 55W 

78 

016 

32.37N 

64.83W 

91 

066 

28.22N 

177.37W 

Sh  i  p 

4yn 

3O.OON 

l4o.OOW 

Ship 

Victor 

34.OON 

164. OOE 

Latitude  35-4S*K 


Latitude  45-S5*N 


Station 

Oegt  lat. 

Deg.  long. 

Block 

Station 

Deg.  lat. 

Deg.  long 

159 

4i.68n 

1.07W 

07 

170 

48.07N 

5.03E 

509 

38.75N 

27.O8W 

07 

354 

46.87N 

1.73E 

308 

36.88N 

76.2OW 

10 

610 

49.95N 

6.57E 

327 

36.12N 

86.68W 

70 

4i4 

52.72N 

I74.IOE 

365 

35.23N 

101 .70W 

70 

454 

5I.88N 

I76.65W 

4o5 

38.85N 

77.03W 

72 

655 

45.58N 

94.I8W 

429 

39.9ON 

84.20W 

72 

712 

46.87N 

68.oew 

445 

38.97N 

92.37W 

72 

764 

46.77N 

100.75W 

451 

37.77N 

99.97W 

72 

775 

47.50N 

111.35W 

476 

39.1  ON 

IO8.53W 

72 

798 

48.80N 

124.73W 

495 

37.73N 

122 .20W 

72 

807 

47.20N 

54.OOW 

506 

4i.25N 

70.07W 

72 

811 

50.22N 

66.27W 

518 

42.75N 

73.8OW 

72 

816 

53.32N 

60.42W 

553 

4i.37n 

96.oew 

72 

836 

51 .27N 

80.65W 

572 

4o.77N 

III.98W 

72 

848 

55.83N 

89.87W 

576 

42.80N 

108. 72W 

72 

879 

53.57N 

113. 52W 

583 

4o.90N 

II7.8OW 

Sh  i  p 

4yc 

52.8ON 

33.5OW 

597 

42.38N 

I22.87W 

662 

44.20N 

IO3.O5W 

4yd 

44. DON 

Latitude  SS-eS’N 


Latitude  65-90“N 


Block 

Stat i on 

Deg.  lat. 

Deg.  long. 

04 

018 

63.98N 

22.63W 

o4 

270 

6I.I8N 

45.42W 

70 

251 

62.97N 

115.62W 

70 

261 

64.82N 

147. 87W 

70 

308 

57.1 5N 

1 70.22W 

70 

350 

57-75N 

152. 52W 

70 

561 

59.52N 

I38.67W 

72 

906 

58.  ION 

68.43W 

72 

907 

58.45N 

78.1 5W 

72 

909 

63.75N 

68.55W 

72 

913 

58.75N 

94.07W 

72 

926 

64.50N 

96.OOW 

72 

954 

60.02N 

111.97W 

72 

945 

58.83N 

122. 58w 

72 

964 

6O.72N 

135.07W 

Ship 

4yb 

56.5ON 

51  .oow 

Block 

Station 

Deg.  lat. 

Deg.  long. 

o4 

202 

76.52N 

68.83W 

70 

026 

71 .30N 

156.78W 

70 

086 

70.12N 

143.67W 

70 

133 

66.87N 

162.63W 

70 

273 

61 .17N 

149.98W 

70 

326 

58.68N 

I56.65W 

72 

917 

80.OON 

85.93W 

72 

924 

74.72N 

94.98W 

72 

938 

67.78N 

115.25W 

72 

968 

68.23N 

135. OOW 

74 

o43 

65.3ON 

126.85W 

74 

051 

72.OON 

124.50W 

74 

072 

76.25N 

119-35W 

74 

074 

78.78N 

103.53W 

74 

082 

82.50N 

62.33W 

74 

090 

70.45N 

68.55W 

APPENDIX  B.  VERTICAL-EXTRAPOLATION  EQUATIONS  * 
B.l  TRC  Equations,  for  500-100  mb 


B.1.1  January 
B.l.  1.1  0-25‘’N 


400  mb 

=  706.59  +  0.93076  (500Zp)  +  : 

400  mb 

^0 

=  -  309.64  +  0.52294  (400Zq)  -  i 

+  0.22918  (400T  ) 

“1^ 

300  mb 

^0 

=  676.43  +  1.6800  (500T^)  +  1 

-  0.86315 (500Z  ) 

JL^ 

300  mb 

^0 

=  984.64  +  0.90022  (400Z^)  +  : 

300  mb 

^0 

=  -  197.35  +  0.527135  (300T  ) 

-  0.275125  (500T  ) 

""  XJi 

300  mb 

^0 

=  -  198.41  +  0.25785  (300Zq)  -  i 

-  0.22293  (400T  ) 

250  mb 

^0 

=  628.86  +  0.95084  (300Z^)  +  : 

250  mb 

^0 

=  -  214.03  +  0.34689 (250Zq)  -  1 

-  0.21443  (300T  ) 

200  mb 

^0 

=  731.60  +  0.95459  (250Z^)  +  : 

200  mb 

'^0 

=  -  405.87  +  0.67632 (200Z^)  -  ( 

150  mb 

^0 

=  1065.2  +  0.91181  (200Z^)  +  : 

150  mb 

^0 

=  -  357.65  +  0.24019  (150T  ■ 

-  0.67665  (200T^) 

100  mb 

^0 

=  1200.9  +  0.85706  (150Z^)  +  ( 

+  1.25420  (150T^) 

100  mb 

^0 

=  -  153.02  +  0.14356 (lOOZ^)  -  1 

B.1.1.2 

25-35“N 

400  mb 

% 

=  631.18  +  0.96996  (500Zq)  +  ; 

400  mb 

% 

=  -  358.96  +  0.60341 (400Z^)  -  ( 

-12' 


0.14444  (500Zq) 


-12' 


-12' 


-12' 


0.46029  (200Zq) 


0.41905 (150Z  ] 


-12' 


0.60341  (500Zq)  -  0.36165  (500T^) 


♦Units  for  all  equations  are  degrees  centigrade  or  tens  of  feet,  as  £q}propriate. 
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300  mb  Zq  =  899.29  +  0.937  58  (5002^^)  +  3.1139  (SOOT^)  +  0.38847  (300Z_j2) 

-  0.38847  (500Z_j^2) 

300  mb  Zp  =  819.14  +  0.96470  (4002^)  +  2.7277  (400T^) 

300  mb  Tq  =  -  339.35  +  0.24933  (3002^^)  -  0.24933  (5002^)  -  0.58345  (5007^^) 

+  0.25733  (300T_j^2) 

300  mb  Tjj  =  -  374.62  +  0.49173  (3002^)  -  0.49173  (4002^^)  -  0.61178  (400Tq) 

+  0.19677  (300T_j^2) 

250  mb  =  592.55  +  0.96311  (3002^)  +  1.9402  (300Tq) 

250  mb  Tjj  =  -  413.03  +  0.85703  (2502^)  -  0.85703  (3002^)  -  0.70181  (300T^) 

+  0.15055 (250T  ) 

*  1a 

200  mb  Zq  =  788.96  +  0.93999  (2502^)  +  2.1985  (250Tq) 

200  mb  Tq  =  -  421.78  +  0.71018  (2002^^)  -  0.71018  (2502^^)  -  0.81701  (250Tjj) 

+  0.17457  (200T_^2) 

150  mb  Zq  =  795.04  +  0.93229  (200Zjj)  +  1.5689  (200Tq)  +  0.26118  (150Z_^2) 

-  0.26118  (200Z_j^2) 

150  mb  =  -  224.03  +  0.36425(1502^)  -  0.38359  (2002^^)  -  0.40238  (200Tq) 

100  mb  Zy  =  1378.7  +  0.85271  (1502^)  +  1.3687  (150Tq)  +  0.23046  (100Z_^2) 

-  0.23046 (1502_^2) 

100  mb  Tq  =  -  147.39  -  0.241283  (1502^)  +  0.21920  (lOOZ^)  -  0.28763  (150Tq) 
+  0.23347 (lOOT  ) 

“1a 

B.1.1.3  35-45’’N 

400  mb  Zp  =  612.98  +  0.97849  (500Zp)  +  2.1177  (500Tq) 

400  mb  Tjj  =  -  356.61  +  0.59834  (4002^)  -  0.59834  (5002^)  -  0.36285  (500Tq) 

300  mb  Zq  =  1036.6  +  0.94251  (500Zjj)  +  3.2431  (500Tq)  +  0.26534  (300Z_^2) 

-  0.26534  (500Z_j^2) 

300  mb  Zq  =  706.65  +  0.95589  (4002^)  +  2.2713  (400Tjj)  +  0.18079  (300Z_j^2) 

-  0.18079  (4002_j2) 

300  mb  Tq  =  -  331.61  +  0.24178  (3002^)  -  0.24178  (5002^^)  -  0.60631  (500?^^) 

+  0.25126  (300T_j^2) 
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300  mb  Tq  =  -  435.68  +  0.57261  (300Z^)  -  0.57261  (4002^^)  -  0.79941  (4007^^) 

+  0.13975  (300T_^2) 

250  mb  =  589.11  +  0.96187 (3002^^)  +  1.8036  (3007^) 

250  mb  Tjj  =  -  520.20  +  1.0741  (2502^^)  -  1.0741  (3002^)  -  1.0116  (300Tjj) 

200  mb  2q  =  653.35  +  0.94888  (2502^^)  +  1.7632  (250Tjj)  +  0.17349  (2002_^2) 

-  0.17349  (2502 

200  mb  T  =  -  235.39  +  0.45738(2002  )  -  0.468919  (2502  )  +  0.23235  (200T  ) 

V  U  V  “  Xa 

-  0.32903  (250Tq) 

150  mb  2  =  572.78  +  0.86206  (2002  )  +  1.3483  (200T  )  +  0.35193(1502  .J 

U  U  V 

-  0.24194(2002  ) 

~*  Xa 

150  mb  Tq  =  -  185.58  +  0.28729  (1502^)  -  0.299338  (2002^)  -  0.22239  (200Tq) 

+  O. 16028 (150T  ) 

**  Xa 

100  mb  2-  =  774.04  +  0.82312(1502.)  +  1.3470  (150T  )  +  0.44376(1002  ) 

V  U  V  *”Xa 

-  0.31899  (1502_^2) 

100  mb  Tq  =  -  317.44  +  0.29877(1002^)  -  0.29877  (1502^)  -  0.40542  (ISOT^) 

+  0.26414  (100T_j^2) 

B.1.1.4  45-55°N 

400  mb  2q  =  583.19  +  0.99329  (5002^)  +  2.0240  (5001^) 

400  mb  Tq  =  -  392.63  +  0.66063(4002^)  -  0.66063  (5002^)  -  0.50001  (SOOT^^) 

300  mb  2  =  983.89  +  0.97663(5002  )  +  2.9589  (500T  )  +  0.24943(3002  ) 

U  U  U  “  Xa 

-  0.24943(5002  ) 

*“  X^ 

300  mb  2  =  698.53  +  0.96027(4002  )+  2.1959  (400T  )  +  0.17041(3002  ) 

U  U  V  ^Xa 

-  0.17041(4002 

300  mb  T  =  -  377.28  +  0.22189(3001  )  +  0.27573(3002  )  -  0.27573(5002  ) 

U  *Xm  u  U 

-  0.78291  (500Tq) 

300  mb  Tp  =  -  459.29  +  0.60299  (3002^)  -  0.60299  (4002^)  -  0.90127  (400T^) 

+  0.13064  (300T_  ^2) 

260  mb  2  =  18.857  +  1.0106(3002  )  +  0.87163(2502  ,„)  -  0.87163(3002  ,.) 

U  U  “  Xa  ••  Xa 
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1 


250  mb 
200  mb  Zq 

200  mb  Tq 

150  mb  Zq 

150  mb  Tq 

100  mb  Zq 

100  mb 


-  5.8706  +  0.20380  (250T  )  +  0.46161  (300T  )  -  0.232416  (300Z  ) 

“  IZ  U  V 

+  0.20080  (250Zjj) 

32.924  +  0.90456  (250Z-)  +  0.83179  (200Z  )  -  0.71065  (250Z  ) 

+  0.70388  (250T^) 

-  9.9353  +  0.30375  (200T_^2)  +  0.30589  (250Tq)  +  0.14406  (200Zq) 

-  0.168146 (250Z^) 

458.71  +  0.91638  (200Z  )  +  1.6339  (200T  )  +  0.34198  (150Z  ) 

U  U  “ 

-  0.25056  (200Z  ) 

-  321.72  +  0.41902  (150Z^)  -  0.41902  (200Zp)  -  0.38313  (200Tjj) 

+  0.12514 (150T  ) 

** 

647.23  +  1.0036  (150ZJ  +  2.1947  (150T.)  +  0.36066  (lOOZ  ) 

0  u  "“IZ 

-  0.36066  (150Z  ) 

**  X^ 

-  295.05  +  0.27568  (lOOZ  )  -  0.27568  (150Z  )  +  0.25378  (lOOT  ) 

U  U  *  Xm 

-  0.34876  (150Tq) 


B.1.1.5 
400  mb  Zq 

400  mb  Tq 

300  mb  Zq 

300  mb  Zq 

300  mb  Tq 

300  mb  Tq 

250  mb  Zq 

250  mb 


so-es’N 


=  602.16  +  0.98297  (500Z^)  +  2.0666  (500T^) 

=  -  373.77  +  0.62642  (400Z^)  -  0.62642  (500Zp)  -  0.4489-1  (SOOT^^) 

=  909.87  +  0.96081  (500Z  )  +  2.6729  (500T  )  +  0.33060  (300Z  ) 

V  V  “  X« 

-  0.33060  (500Z  ) 

~  Xm 

=  625.64  +  0.96621  (400Z  )  +  1.9379  (400T  )  +  0.24842  (300Z  ) 

V  U  ”*X^ 

-  0.24842  (400Z  ) 

=  -  310.23  +  0.36138  (300T  )  +  0.22666  (300Z  )  -  0.22666  {500Z  ) 

*  Xm  V  V 

-  0.63937  (500Tq) 

=  -  444.93  +  0.58513  (300Zp)  -  0.58513  (400Zq)  -  0.88524  (400Tjj) 

+  0.18519  (300T  ) 

~  Xm 

=  447.70  +  0.97127  (300Z  )  +  1.5279  (300T  )  +  0.25606  (250Z  ) 

V  V  “  X^ 

-  0.25606  (300Z  ) 

“*1^ 

=  -  503.12  +  1.0475  (250Zp)  -  1.0475  (300Zjj)  -  0.87663  (300T^) 
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200  mb  Zq  =  375.34  +  0.91231  (2502^)  +  1.4037  (250Tq)  +  0.399646(2002 

-  0.31708(2502  ) 

* 

200  mb  Tq  =  -  363.01  +  0.61512  (2002^)  -  0.61512  (2502^)  -  0.52050  (250Tq) 

+  0.14714  (200T  ) 

150  mb  2^  =  407.21  +  0.93279(2002^)  +  1.6998  (200Tq)  +  0.334942  (150Z_^2) 

-  0.24387(2002  ) 

“1^ 

150  mb  T  =  -  209.58  +  0.27673(1502  )  -  0.27673  (2002  )  +  0.425375  (150T  ) 

U  U  U  “  Xa 

-  0.199655  (200T_j^2) 

100  mb  2^  =  649.40  +  0.3378293(1002  +  1.0090  (1502^)  +  2.5856  (150T^) 

-  0.33435(1502  ) 

“  Xtt 

100  mb  T  =  -  149.51  +  0.14097(1002 J  -  0.14097  (1502  )  +  0.47163(100T  ) 

0  U  1/  “  x^ 


B.1.1.6 
400  mb  2^ 
400  mb 
300  mb  2^ 

300  mb  2^ 

300  mb 

300  mb  Tq 

250  mb  Zq 

250  mb  Tq 
200  mb  2^ 

200  mb  Tq 


65-90°N 

605.56  +  0.98132(5002^)  +  2.0803  (500Tq) 

-  389.46  +  0.65392  (4002 jj)  -  0.65392  (5002^^)  -  0.50034  (SOOT^) 
867.97  +  0.96473  (5002  J  +  2.5628  (500T  )  +  0.35792(3002 

-  0.35792(5002  ,„) 

635.27  +  0.96765  (4002^)  +  1.9720  (400T^)  +  0.22945  (300Z_^2) 

-  0.22945(4002  ) 

““X^ 

-  338.38  +  0.33706  (300T  )  +  0.24  892(3002  )  -  0.24892(5002  ) 

Xu  V  V 

-  0.68280  (500T^) 

-  452.77  +  0.59744  (300Zp)  -  0.59744  (4002^^)  -  0.87503  (400Tq) 

+  0.17470  (300T 

414.74  +  0.97438  (300Zq)  +  1.4583  (SOOT^)  +  0.31010  (250Z_^2) 

-  0.31010(3002 

-  514.61  +  1.0745  (250Zjj)  -  1.0745  (3002^)  -  0.89574  (300Tjj) 

431.24  +  a99334  (2502.)  +  1.7106 (250TJ  +  0.32688(2002  ) 

0  V  ““  Xu 

-  0.32688(2502  ) 

~  Xu 

-  368.27  +  0.62636  (200Zq)  -  0.62636  (2502^)  -  0.57128  (250Tq) 

+  0.19900  (200T 
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150  mb 


=  421.31  +  1.0239  (200Zjj)  +  1.8120  (200Tq)  +  0.32143  (150Z_j2) 

-  0.32143  (200Z  ) 

150  mb  =  -  392.57  +  0.47331  (150Zq)  -  0.46425  (200Zq)  -  0.33403  (200Tq) 

100  mb  Zq  =  558.31  -  0.446846  (100Z_^2)  1.0444  (150Zq)  +  2.3500  (150Tq) 

-  0.42155  (150Z  ) 

100  mb  Tq  =  -  447.15  +  0.42088  (lOOZ^^)  -  0.42088  (WOZ^^)  -  0.61475  (ISOT^) 

B.1.2  April 

B.1.2.1  0-25‘‘N 

400  mb  Z^  =  815.63  +  0.87563  (5002^)  +  2.4202  (500Tq) 

400  mb  =  -  285.00  -  0.48136(5002^)  -  0.24025  (500T^)  +  0.23962  (400T_j^2) 

+  0.48136  (400Zq) 

300  mb  Z„  =  883.28  +  2.9994 (500TJ  +  0.90680(5002.)  +  0.45035(3002  ,„) 

0  0  u 

-  0.45035(5002  ) 

*  xz 

300  mb  Zq  =  1120.8  +  0.84628  (4002^^)  +  3.0236  (4001^) 

300  mb  Tq  =  -  219.00  -  0.15830(5002^)  +  0.392115  (300T_^2)  ~  0.246665  (500T_^2) 

+  0.15830(3002^) 

300  mb  T„  =  -  187.83  -  0.23942  (400Z„)  +  0.28389  (300T  )  +  0.23942  (3002  ) 

250  mb  2^  =  893.50  +  0.87126(3002^)  +  2.2120  (SOOT^^) 

250  mb  =  -  96.865  -  0.18348  (3002^)  +  0.33633  (300Tjj)  +  0.241864  (250T_j^2) 

.  +  0.18348  (2502^) 

200  mb  Zp  =  906.75  +  0.91056(2502^)  +  2.4416  (2507^^) 

200  mb  Tq  =  -  363.48  -  0.60307(2502^)  -  0.43571  (2507^^)  +  0.60307  (2002^) 

150  mb  Zq  =  1218.3  +  0.88270(2002^)  +  2.8453  (2001^) 

150  mb  Tq  =  -  455.60  -  0.58532  (2002^)  -  0.82914  (200T^)  +  0.58532  (1502^) 

100  mb  2^  =  1672.5  +  0.78825  (1502 J  +  1.4079  (150TJ  +  0.23681(1002  ) 

0  u  u  “Xtt 

-  0.23681  (150Z_^2) 

100  mb  Tq  =  -  2,2489  -  0.08623  (1502^^)  +  0.40531  (100T_^2)  +  0.11484  (1002^^) 
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B.1.2.2 

25-35‘N 

400  mb 

^0 

=  632.86  +  0.96998  (500Zjj)  + 

400  mb 

^0 

=  -  240.91  -  0.40183  (500Z^)  + 

300  mb 

=  960.60  +  0.96139  (500Z^)  + 

-  0.30711  (500Z_^2) 

300  mb 

=  809.52  +  0.96985  (400Z^)  + 

300  mb 

^0 

=  -  309.87  -  0.22492  (500Z^)  - 

+  0.22492 (300Zq) 

300  mb 

^0 

=  -  401.92  -  0.52322  (400Zq)  - 

250  mb 

=  596.65  +  0.96130  (300Zp)  + 

250  mb 

^0 

=  -  433.55  -  0.896  89  (SOOZjj)  - 

200  mb 

^0 

=  820.61  +  0.93291  (250Z^)  + 

200  mb 

^0 

=  -  525.65  -  0.88733  (250Zq)  - 

150  mb 

=  926.05  +  0.90937  (200Zq)  + 

-  0.21785 (200Z  ) 

*  Xm 

150  mb 

^0 

=  -  239.56  -  0.393638  (200Zp)  ■ 

+  0.37750  (150Z^) 

100  mb 

% 

=  739.05  +  0.81544  (150Zp)  + 

-  0.38861  (150Z  ) 

~  X^ 

100  mb 

T 

0 

=  -  214.98  -  0.20207  (150Zp)  - 

+  0.20207  (lOOZp) 

B.l. 

,2.3 

35-45'’N 

400  mb 

^0 

=  616.57  +  0.97733  (500Zp)  + 

400  mb 

^0 

=  -  231.14  -  0.38392  (500Zp)  + 

300  mb 

^0 

=  894.67  +  0.99456  (500Zp)  + 

-  0.30064  (500Z  ) 

—  Xfc 

300  mb 

"0 

=  825.63  +  0.95968  (400Zq)  + 

2.3071  (500Tq) 

0.13241  (400T  )  +  0.40188  (400Z-) 

u 

3.6596 (500T)  +  0.30711 (300Z  ) 

U 

2.8255  (400T) 
u 

0.37600  (500T  )  +  0.18817 (SOOT  ) 

0.54745  (400Tp)  +  0.52322  (SOOZ^) 
1.8999  (300Tjj) 

0.62931  (300T^)  +  0.89689  (250Zjj) 
2.3373  (250Tjj) 

1.0234  (250Tq)  +  0.88733  (200Z^) 
1,8297  (200Tq)  +  0.21785  (150Z_j2) 

0.47937  (200Tjj)  +  0.13365  (150T_  12) 

1.0642  (150T^)  +  0.51238  (100Z_j2) 

0.20973  (150T  )  +  0.43163  (lOOT  ) 

V  **  Xa 


2.1962  (500Tp) 

0.13851  {400T  )  +  0.38392  (400Z  ) 

^  x«  u 

3.0219  (500T  )  +  0.30064  (300Z  ) 

U  “  Xtt 

2.5826  (400T^) 
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300  mb  Tq  =  -  399.48  -  0.29395  (SOOZ^)  -  0.81287  (SOOT^)  +  0.19631  (300T_^2) 
+  0.29395  (300Zj|) 

300  mb  Tq  =  -  371.55  -  0.48501  (400Zq)  -  0.62055  (400Tq)  +  0.17766  (300T_^2) 
+  0.48501  (300Zq) 

250  mb  Zq  =  612.07  +  0.95534  (300Zq)  +  1.8600  (300Tjj) 

250  mb  Tq  =  -  528.14  -  1.0919  (OOOZ^)  -  1.0318  (SOOT^)  +  1.0919  (250Zq) 

200  mb  Zq  =  867.73  +  0.91451  (250Z^j)  +  2.0612  (250T^) 

200  mb  Tq  =  -  525.97  -  0.88448  (250Zq)  -  1.0571  (250T^)  +  0.88448  (200Zq) 

150  mb  Zq  =  916.55  +  0.91271  (200Zq)  +  1.6037  (200Tq)  +  0.18821  (150Z_^2) 

-  0.18821  (200Z_j^2) 

150  mb  Tq  =  -  223.62  -  0.37343  (200Z^)  -  0.32115  (200Tq)  +  0.35678  (ISOZ^) 

100  mb  Zq  =  424.07  +  0.78057  (150Z^)  +  1.0429  (ISOTjj)  +  0.61291  (100Z_^2) 

-  0.40261  (150Z_j^2) 

100  mb  Tq  =  -  223.39  -  0.20925  (ISOZ^^)  -  0.19717  (ISOT^)  +  0.37145  (100T_^2) 
+  0.20925  (IOOZq) 


B.  1.2.4  45-55°N 

400  mb  Zq  =  595.19  +  0.98740(5002^)  +  2.0838  (500Tq) 

400  mb  =  -  350.91  -  0.58825(5002^)  -  0.34646  (SOOT^)  +  0.58825(4002^) 

300  mb  Zq  =  993.16  +  0.96966  (500Zjj)  +  3.0018  (500Tp)  +  0.25408  (300Z_j^2) 

-  0.25408  (500Z_^2) 

300  mb  Zq  =  848.13  +  0.94668(4002^)  +  2.4260  (400T^) 

300  mb  Tq  =  -  354.80  +  0.25815(3002^)  -  0.25815  (5002^)  -  0.79333  (500T^j) 

+  0.25774  (300T  „) 

300  mb  Tp  =  -  427.46  -  0.55874(4002^)  -  0.86039  (400Tq)  +  0.17012  (300T_j^2) 
+  0.55874(3002^) 

250  mb  Zq  =  655.12  +  0.93892(3002^)  +  1.7682  (300T^) 

250  mb  =  -  525.48  +  1.0845  (2502^)  -  1.0845(3002^)  -  1.0390  (300T^) 

200  mb  2q  =  854.43  +  0.91151  (2502^^)  +  1.6549  (2501^) 
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200  mb  =  -  464.82  +  0.78173  (2002^^)  -  0.78173  (2502^^)  -  0.81851  (250Tq) 

150  mb  =  896.52  +  0.95176(2002^)  +  1.9208  (200Tp) 

150  mb  Tq  =  -  281.18  -  0.36303(2002^)  -  0.26345  (200Tq)  +  0.12394  (150T  ^g) 

+  0.36303(1502^) 

100  mb  2p  =  885.48  +  0.84649  (1502 p)  +  1.5778(1507^)  +  0.365949(1002 

-  0.27128(1502  ) 

~  XiW 

100  mb  =  -  199.14  -  0. 18391  (1502^)  +  0.41103(1007  "  0.19316(1507 

+  0.18391(1002^) 


B.  1.2.5  55-65°N 

400  mb  2^  =  598.13  +  0.98638(5002^)  +  2.1288(5007^) 

400  mb  =  -  257.67  -  0.42527  (5002^)  +  0.42527  (4002^) 

300  mb  Z  =  906.42  +  0.96076  (5002  )  +  2.7817  (5007  )  +  0.33693(3002  ) 

U  00  “12 

-  0.33693  (500Z  ) 

“1^ 

300  mb  2^  =  668.67  +  0.94384  (4002^)  +  2.0223(4007^)  +  0.26720(3002  ^g) 

-  0.26720(4002  ) 

300  mb  =  -  339.06  +  0.24697(3002^)  -  0.24697  (5002^)  -  0.79880  (5007^) 

+  0.33200(3007  ) 

300  mb  Tp  =  -  480.17  -  0.63012(4002^)  -  1.0139  (4007^)  +  0.13625(3007 
+  0.63012(3002^) 

250  mb  Zq  =  642.20  +  0.94091  (300Zp)  +  1.6371  (3007p) 

250  mb  =  -  480.82  -  0.98927  (3002^^)  -  0.89260  (3007^)  +  0.98927  (2502^) 

200  mb  2^  ■-=  743.49  +  0.94421  (2502^)  +  1.5948  (2507^) 

200  mb  Tq  =  -  395.31  -  0.66475  (2502 p)  -  0.55102  (2507^)  +  0.66475(2002^) 

150  mb  Zq  =  799.77  +  0.98016  (200Zjj)  +  2.1219  (2007^) 

150  mb  Tjj  =  -  211.34  +  0.27383  (2002^)  +  0.305215(1507  ^2)  “  0.156225(2007 

+  0.27383  (150Zq) 

100  mb  2^  =  963.81  +  0.94337  (1502^^)  +  1.8865  (1507^)  +  0.28279(1002 

-  0.28279(1502  ) 

“X^ 

100  mb  '’^0  =  ~  71.404  -  0.1064059  (1502^)  +  0.40428(1007  +  0.096748  (1002^) 
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B.1.2.6  eS-QO^N 


400  mb 

^0  = 

609.41  +  0.97852  (500Zq)  +  2.0590  (SOOT^) 

400  mb 

II 

O 

H 

-  471.71  -  0.79801  (500Z^)  -  0.79563  (500Tq)  +  0.79801  (400Zjj) 

300  mb 

725.98  +  0.94652  (500Zq)  +  1.9040  (5007^)  +  0.49395  (300Z_^2) 

-  0.49395  (500Z_j2) 

300  mb 

626.44  +  0.93500  (400Zq)  +  1.8018  (400Tjj)  +  XL35257  (300Z_j^2) 

-  0.35257  (400Z  ) 

^  IZ 

300  mb 

^0  = 

-  326.71  -  0.2311103  (500Zp)  -  0.81894  (5007^)  +  0.34474  (300T_j^2) 

+  0.23308  (300Zp) 

300  mb 

-  428.44  -  0.55586  (400Zq)  -  0.92731  (400Tq)  +  0.17313  (300T_j2) 

+  0.55586  (300Zq) 

250  mb 

N 

o 

II 

596.87  +  0.95509  (300Zq)  +  1.5431  (300Tq) 

250  mb 

^0  = 

-  464.32  -  0.95773  (300Zp)  -  0.80865  (3007^)  +  0.95773  (250Zq) 

200  mb 

^0  = 

662.76  +  0.97083  (250Z^)  +  1.7246  (2507^) 

200  mb 

^0  = 

-  376.28  -  0.63434  (250Zq)  -  0.46781  (2507^)  +  0.63434 (ZOOZ^j) 

150  mb 

814.12  +  0.97731  (200Zjj)  +  2.1865(2007^) 

150  mb 

o 

II 

-  318.99  -  0.42124  (200Zjj)  -  0.34033(2007^)  +  0.15570  (150T_j2) 

+  0.42124  (150Z^) 

100  mb 

II 

o 

N 

728.40  +  0.95857  (150Zq)  +  1.7110  (1507^^)  +  0.46780  (100Z_^2) 

-  0.46780 (150Z_^2) 

100  mb 

^0  = 

-  79.792  -  0.1142054  (150Zq)  +  0.44908  (1007  ^^)  +  0.10534  (lOOZ^^) 

B.1.3  Jiily 
B.  1.3.1  0-25*N 

300  mb  Zq  =  1869.7  +  0.80469  (SOOZ^)  +  3.0136  (SOOT^^)  +  0.17441  (300Z_j2) 

-  0.40373  (500Z_j^2) 

300  mb  T„  =  -  244.51  -  0.17691  (500ZJ  +  0.307235  (300T  -  0.166365  (SOOT  ) 

0  U  •’Itt  “AZ 

+  0.17691  (300Zjj) 

=  1021.7  +  0.82945  (300Z^)  +  2.0387  (3001^^) 
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250  mb 


\ 


250  mb  T.  «  -  57.322  -  0.10317(3002.)  +  0.64580 (300TJ  +  0.16236 (250T 
0  U  U 

+  0.10317  (250Zjj) 

200  mb  Zq  =  1221.5  +  0.82421  (250Zjj)  +  2.5376  (250Tq) 

200  mb  Tq  =  -  251.68  -  0.40912  (2502^)  +  0.40912  (2002^) 

150  mb  Zq  =  1529.3  +  0. 80811  (2002^^)  +  2.9565  (200T^) 

150  mb  Tjj  =  -  665.75  -  0.63957  (2002^^)  -  1.2876  (200Tq)  +  0.67250  (1502^) 

100  mb  2.  =  933.76  +  0.85394(1502.)  +  0.67683(1002  ,.)  -  0.67683(1502  ,,) 

100  mb  Tp  =  -  119.88  -  0.10633(1502^)  +  0.50032  (100T_j^2)  +  0.10633  (1002^) 

B.  1.3.2  25-35'’N 

300  mb  2.  =  726.53  +  1.0038  (5002  J  +  2.3720  (500T.)  +  0.42381(3002  ) 

0  0  U 

-  0.42381  (500Z_j^2) 

300  mb  T„  =  -  208.38  -  0.14939(5002.)  +  0.371545  (300T  -  0.225995  (500T  ,, 

+  0.14939 (300Zq) 

250  mb  Zq  =  781.63  +  0.90467  (3002^^)  +  1.9867  (300T^) 

250  mb  =  -  150.38  -  0.29403  (3002^^)  +  0.31531  (250T_^2)  0.29403  (2502^^) 

200  mb  2^  =  889.39  +  0.91225  (2502^)  +  2.1325  (250Tjj) 

200  mb  Tq  =  -  423.73  -  0.70808(2502^)  -  0.66739  (2501^)  +  0.70808(2002^^) 

150  mb  Zjj  =  1337.5  =  0.85442  (200Zjj)  +  2.8475 (200Tq) 

150  mb  Tq  =  -  462.85  -  0.59311  (2002^)  -  0.89370  (2001^)  +  0.59311  (1502^) 

100  mb  2^  =  497.76  +  0.94345  (1502^)  +  0.70659  (100Z_^2)  -  0.70659  (150Z_^2) 

100  mb  T„  =  -  169.13  -  0.14842  (1502 J  -  0.25307  (150T.)  +  0.50227  (lOOT 

0  0  0 

+  0.14842  (lOOZjj) 

B.  1.3.3  35-45°N 

300  mb  Zq  =  847.16  +  0.96789  (5002^)  +  3.4309  (500T^)  +  0.38817  (300Z_^2) 

-  0.38817  (500Z_j^2) 

300  mb  Tq  =  -  372.81  -  0.27089  (5002^)  -  0.35466  (500T^)  +  0.27089(3002^) 

250  mb  Zq  =  622.67  +  0.95408  (300Zq)  +  1.9574  (300Tq) 

250  mb  Tp  =  -  208.16  -  0.42348  (3002^)  +  0.21194  (250T_j^2)  +  0.42348(2502^) 
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200  mb  =  849.64  +  0.92393(2502^)  +  2.2085  (2501^) 

200  mb  Tq  =  -  494.01  -  0. 82225 (2502^^)  -  1.0240  (250Tq)  +  0.82225  (2002^^) 

150  mb  2q  =  1407.9  +  0.82718  (2002 +  2.0969  (200Tjj) 

150  mb  Tq  =  -  355.14  -  0.45799(2002^)  -  0.65922  (200Tq)  +  0.23149  (150T_^2) 

+  0.45799(1502^) 

100  mb  2q  =  492.97  +  0.76505  (1502^^)  +  1.3464  (150Tq)  +  0.61936  (1002_^2) 

-  0.40620(1502  ) 

100  mb  Tq  =  -  134.84  -  0.12303(1502^)  +  0.46319  (100T_j^2)  0.12303  (1002 p) 

B.1.3.4  45-55°N 

300  mb  2^  =  1316.5  +  0.97715(5002^)  +  4.7636  (500T^) 

300  mb  =  -  425.59  -  0.31064  (5002^)  -  0.69662  (500Tp)  +  0.31064  (3002^) 

250  mb  2q  =  616.80  +  0.95413  (3002 p)  +  1.8564  (SOOT^) 

250  mb  Tq  =  -  514.49  -  1.0602  (3002^^)  -  1.0252  (300T|j)  +  1.0602  (2502^) 

200  mb  2^  =  914.67  +  0.90037  (250Zp)  +  1.8947  (250T^) 

200  mb  =  -  345.11  -  0.646616  (2502^)  -  0.72702  (250T^)  +  0.11790  (200T_j^2) 
+  0.63537  (2002^) 

150  mb  2^  =  1185.9  +  0.87478  (2002^)  +  1.5345  (200T^) 

150  mb  Tjj  =  -  194.32  -  0.375396  (2002^)  -  0.34041  (200Tp)  +  0.35238  (1502^) 

100  mb  Z„  =  941.26  +  0.80547(1502.)  +  1.5611  (150TJ  +  0.37735(1002  ) 

0  U  V  ”XZ 

-  0.25469(1502 

100  mb  =  -  163.50  -  0.15029  (1502^)  +  0.34838  (100T_j^2)  +  0.15029  (1002^) 

B.1.3.5  55-65°N 

300  mb  2^  =  1303.3  +  0.98361  (5002^)  +  4.6744  (500T^) 

300  mb  T„  =  -  342.10  -  0.25104  (5002  )  -  0.69667  (500T  )  -  0.29083  (300T  ) 

0  U  u  ~x^ 

+  0.25104  (300Zq) 

250  mb  2p  =  667.66  +  0.93693  (3002^^)  +  1.8452  (300T^) 

250  mb  Tp  =  -  484.25  -  1.0001  (3002^)  -  1.0449  (300T^)  +  0.14462  (250T_j^2) 

+  1.0001  (2502^) 
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200  mb  Zq  =  991.66  +  0.87344 (250Zjj)  +  1.6035  (250Tq) 

200  mb  Tq  =  -  374.10  -  0.62516  (250Zjj)  -  0.66289  (250Tq) 

+  0.62516  (200Z^) 

150  mb  Zq  =  869.92  4  0.95889  (200Z^)  +  1.7584  (200Tq) 

150  mb  Tq  =  -  294.43  -  0.37650  (200Z^)  -  0.25674  (200T^) 

100  mb  Zq  =  733.90  +  0.89688  (150Zjj)  +  2.7989  (150T^)  H 

100  mb  Tq  =  -  93.496  -  0.Q80453  (1502^^)  +  0.35023  (150Tq: 

-  0.207775  (150T  )  +  0.080453  (lOOZ  ) 

JLa  u 


B.  1.3.6 
300  mb  Zq 

300  mb  Tq 

250  mb  Zq 
250  mb  Tq 

200  mb  Zq 
200  mb  Tq 

150  mb  Z^ 
150  mb 
100  mb  Zq 
100  mb  Tq 


65-90°N 

1110.4  +  0.94325  (500Zjj)  +  3.8614  (500T^)  H 

-  0.21206  (500Z  ) 

-  398.74  -  0.29542(5002^)  -  0.987281  (SOOT^; 
+  0.29542  (300Zq) 

703.77  +  0.92276(3002^)  +  1.7024  (SOOT^^) 

-  451.46  -  0.92596  (300Zp)  -  0.96343  (300T^) 
+  0.92596  (250Zq) 

1008.9  +  0.86339(2502^)  +  1.2375  (2501^) 

-  310.26  -  0.51682  (2502^)  -  0.46979  (250Tp) 
+  0.51682(2002^) 

850.21  +  0.96417  (2002^)  +  1.7373  (200T^) 

-  220.80  -  0.27929  (200Zp)  +  0.27929  (1502^) 
561.38  +  0.91608(1502^)  +  3.0046  (150T^)  H 

-  72.710  -  0.066889  (150Zp)  +  0.40139  (150T^; 

-  0.21265  (150T  )  +  0.066889  (1002  ) 

"  Aa  V 


B.1.4  October 
B.  1.4.1  0-25°N 

400  mb  Zq  =  749.46  +  0.90882  (5002 ij)  +  1.9588  (500Tq) 

400  mb  Tq  =  -  267.28  -  0.45034  (5002^^)  -  0.22181.(500Tjj) 

+  0.45034  (400Zq) 


+  0.14864  (ZOOTjg) 

+  0.37650  (150Zq) 

-  0.13696  (lOOZ  j^g) 

I  +  0.375615  (lOOT  ) 

**  Am 

■  0.21206  (300Z_j^2) 

I  +  0.32780  (300T  ) 

^  Am 

+  0.14414 (250T  ) 

“Am 

+  0.18578 (200T  ) 

"  Am 


0.15499(1002  ) 

~  Xm 

+  0.51096  (lOOT  ) 

“Am 


+  0.24780  (400T  ) 

“  Am 
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300  mb 

"o 

=  922.57  +  0.89896(5002^)  + 

-  0.43019(5002 

300  mb 

% 

=  1199.5  +  0.81466 (4002q)  + 

300  mb 

^0 

=  -  253.87  -  0.18387  (5002^)  - 

+  0.18387(3002^) 

300  mb 

^0 

=  -  198.82  -  0.25477  (4002jj)  + 

250  mb 

=  980.13  +  0.84434(3002^)  + 

250  mb 

^0 

=  -  77.465  +  0.14500  (2502 jj)  - 

+  0.20815  (2507_^2) 

200  mb 

^0 

=  1183.6  +  0.83799  (2502^)  + 

200  mb 

^0 

=  -  235.01  -  0.38991(2502^)  + 

+  0.38991  (2002^) 

150  mb 

^0 

=  1445.7  +  0.83093  (2002^)  + 

150  mb 

^0 

=  -  413.31  -  0.54003  (2002jj)  - 

+  0.54003  (1502^) 

100  mb 

^0 

=  1200.4  +  0.86289  (1502^)  + 

-  0.40167  (1502_j2) 

100  mb 

^0 

=  -  104.45  +  0.096897  (1002^)  - 

B.  1.4.2 

25-3  S'N 

400  mb 

^0 

=  621.82  +  0.97594  (5002p)  + 

400  mb 

^0 

=  -  386.06  -  0.65129(5002^)  - 

300  mb 

=  817.16  +  0.94761  (500Zp)  + 

-  0.44116(5002  ) 

^  Xa 

300  mb 

=  877.32  +  0.94372  (4002p)  + 

300  mb 

^0 

=  -  272.55  -  0.19914  (SOOZjj)  - 

+  0.19914  (300Zq) 

300  mb 

^0 

=  -  288.62  -  0.37651  (400Zjj)  - 

+  0.37651  (3002q) 

2.6759  (500T  )  +  0.43019(3002 

V 

2.9930  (400Tq) 

0.21738  (500T)  +  0.29644  (300T  ) 

U  •  Xtt 

0.26961(300T  +  0.25477  (3002  J 

2.2405  (300Tjj) 

0.14500  (3002  +  0.49616  (300Tj^) 

2.8257  (250Tjj) 

0.303905(2007  -  0.216295(2507 

3.1419  (2007jj) 

0.73712(2007  )  +  0.15879(1507  ) 

u  *•*<« 

1.4277(1507^)  +  0.40167  (1002_j^2) 
0.096897(1502  )  +  0.61744(1007  ) 

V  “  Xm 

2.2969  (5007^) 

0.44416(5007^)  +  0.65129  (4002^) 
3.0894(5007  )  +  0.44116(3002  ) 

V  *~Xm 

2.9377  (4007^) 

0.38498(5007  )  +  0.33833(3007  ) 

U  “  Xa 

0.32711(4007.)  -  0.25106(3007 

V  “•ix 
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250  mb  Zq  =  596.27  +  0.96151  (3002^)  +  1.8806  (SOOT^) 

250  mb  =  -  220.45  -  0.45167(3002^)  +  0.33843  (2 50T_^2)  -  0.18064  (300T_^2) 
+  0.45167  (2502^) 

200  mb  Zjj  =  766.32  +  0.94682  (250Zjj)  +  2.1916  (250Tq) 

200  mb  T  =  -  325.87  -  0.54177(2502  )  -  0.55888 (250T  )  +  0.24364(200T  ) 

V  U  U 

+  0.54177  (200Zq) 

150  mb  Zq  =  1091.0  +  0.90587  (2002 jj)  +  2.2040  (200Tjj) 

150  mb  T  =  -  393.31  -  0.51311(2002  )  -  0.71958 (200T  )  +  0.20260(150T  ) 

U  U  U  ■*  IZ 

+  0.51311  (1502q) 

100  mb  2  =  1184.8  +  0.87474(1502  J  +  1.6587  (150TJ  +  0.37570(1002  ) 

U  U  U 

-  0.37570(1502  ) 

100  mb  T  =  -  149.61  -  0.14316(1502  )  +  0.49699  (lOOT  )  +  0.14316(1002  ) 

U  U  ^Xi^  u 

B.  1.4.3  35-45‘N 

400  mb  2q  =  632.57  +  0.96932  (5002^)  +  2.2196  (500Tp) 

400  mb  =  -  344.03  -  0.57721  (5002^)  -  0,30795  (500Tp)  +  0.57721(4002^) 

300  mb  2^  =  1008,2  +  0.95621  (500Zjj)  +  3.4029  (500Tq)  +  0.27219  (300Z_j^2) 

-  0.27219  (500Z_j2) 

300  mb  2^  =  808,53  +  0.96890  (4002jj)  +  2.7007  (400Tq) 

300  mb  T  =  -  355,77  -  0.26194(5002  )  -  0.59606  (500T  )  +  0.23115  (300T  ) 

U  U  U  “  Xm 

+  0.26194  (3002q) 

300  mb  T  =  -  395.58  -  0.52021(4002  )  -  0.66210  (400T  )  +  0.17812  (300T  ) 

V  V  V  *Xa 

+  0.52021  (300Zq) 

250  mb  2^  =  605.70  +  0.95904  (3002^)  +  1.9353  (300Tp) 

250  mb  T  =  448,66  -  0.93017(3002  )  =  0.82981  (300T  )  +  0.12364(250T  ,„) 

+  0.93017  (250Zq) 

200  mb  2^  =  767.94  +  0.94484  (2502^)  +  2.0803  (250Tjj) 

200  mb  Tjj  =  -  479.16  -  0.80121  (2502^)  -  0.90957  (250T^)  +  0.80121  (2002^) 

150  mb  2^  =  1041.6  +  0.91377(2002^)  +  1,8716  (200Tq) 

150  mb  Tjj  =  -  285.18  -  0.446783(2002^)  -  0.48940  (2001^)  +  0.43194  (1502^) 
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100  mb 


=■  656.20  +  0.82111  (150Z^)  +  1.4108  (ISOT^)  +  0.53913  (lOOZ 

-  0.40248  (150Z  ) 


100  mb 

^0  = 

-  252.24  -  0.23388  (150Z  )  -  0.29674  (150T  )  +  0.33362  (lOOT  ) 

u  U  "  Xa 

+  0.23388  (IOOZq) 

B.1.4.4 

45-55'’N 

400  mb 

=^0  = 

595.32  +  0.98839  (500Z^)  +  2.1288  (500Tq) 

400  mb 

^0  = 

-  256.60  -  0.42436  (500Zjj)  +  0.42436  (400Zp) 

300  mb 

1090.4  +  0.96440  (500Z^)  +  3.5450  (500T^)  +  0.19219  (300Z 
-  0.19219 (500Z  ) 

300  mb 

^0  = 

844.40  +  0.95304(4002^)  +  2.6843  (400Tq) 

300  mb 

^0  = 

-  376.03  -  0.27651  (500Zq)  -  0.73220  (500Tjj)  +  0.22616  (SOOT 
+  0.27651  (300Zq) 

300  mb 

^0  = 

-  413.10  -  0.54190(4002^)  -  0.73773  (400Tp)  +  0.15365(3001 
+  0.54190  (300Zq) 

250  mb 

596.22  +  0.95972(3002^)  +  1.8002  (300T^) 

250  mb 

^0  = 

-  507.64  -  1.0474(3002^)  -  0.97244  (SOOT^)  +  1.0474  (.2502^) 

200  mb 

^0  = 

768.52  +  0.94137(2502^)  +  1.8987  (250T^) 

200  mb 

^0  = 

-  25.041  -  0.202209(2502^)  +  0.25833  (250Tp)  +  0.29553  (200T 
+  0.17758(2002^) 

150  mb 

987.15  +  0.92914(2002^)  +  2.0046  (200T^) 

150  mb 

^0  = 

-  86.360  -  0.256806  (2002p)  +  0.22912  (1502 j^) 

100  mb 

773.43  +  0.82019(1502^)  +  1.5487  (150Tq)  +  0.43879(1002 

-  0.30660  (1502_j^2) 

100  mb 

-  346.47  -  0.32148  (1502p)  -  0.32855  (150T^)  +  0.32148(1002^) 

1.4.5 

55-65°N 

400  mb 

584.54  +  0.99394  (5002p)  +  2.1315  (SOOT^) 

400  mb 

^0  = 

-  362.53  -  0.61015  (5002p)  -  0.35196  (500T^)  +  0.61015  (4002^) 

300  mb 

^0  = 

1023.6  +  0.98663  (5002jj)  +  3.4503  (SOOT^)  +  0.21071(3002 
-  0,21071(5002  ^g) 
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300  mb  =  797.15  +  0.97012  (4002^)  +  2.5359  (400Tjj) 

300  mb  Tg  =  -  369.11  -  0.27178  (5002^)  -  0.78274  (500T^j)  +  0.28637  (300T_j2) 

+  0.27178  (300Zq) 

300  mb  Tq  =  -  410.73  -  0.53915  (4002^)  -  0.77094  (400T^)  +  0.19345  (300T_j2) 
+  0.53915  (300Zq) 

250  mb  Zq  =  618.35  +  0.95092  (3002^)  +  1.7531  (300Tq) 

250  mb  Tq  =  -  496.94  -  1.0262  (3002^)  -  0.92485  (300Tq)  +  1.0262  (2502^) 

200  mb  Z^  =  763.87  +  0.93930  (2502^)  +  1.7082  (250Tq) 

200  mb  =  430.17  -  0.72426(2502^)  -  0.67308  (250T^)  +  0,72426  (2002^) 

150  mb  2^  =  864.45  +  0.96060  (2002^)  +  1.9861  (200Tq) 

150  mb  =  -  346.83  -  0.45282  (2002^)  -  0,34739  (200Tq)  +  0.45282  (1502^) 

100  mb  2q  =  742.88  +  0.98276  (1502^)  +  2.1891  (150Tq)  +  0.35666  (1002_^2) 

-  0.35666 (1502_^2) 

100  mb  Tq  =  -  204.30  -  0.19001  (1502^)  +  0.22212  (100T_j^2)  +  0.19001(1002^) 
B.  1.4.6  65-90°N 

400  mb  Zq  =  622,47  +  0.97229  (6002^^)  +  2.1412  (500Tq) 

400  mb  Tq  =  -  346.70  -  0.58108  (5002^)  -  0.31663  (500?^)  +  0.58108  (4002^) 

300  mb  Z^  =  1079.0  +  0.92800  (5002^)  +  3.2584  (500T^)  +  0.24573  (300Z_^2) 

-  0.24573(5002  ) 

300  mb  Zq  =  861.41  +  0.94166  (4002^^)  +  2.5871  (400Tq) 

300  mb  Tq  =  -  410.43  -  0.30303  (5002^)  -  0.88779  (500Tp)  +  0.24852  (300T_^2) 

+  0.30303  (3002^) 

300  mb  Tq  =  -  535.63  -  0.70453  (4002^^)  -  1.0405  (400Tq)  +  0.70453  (3002^) 

250  mb  Zq  =  632.52  +  0.94520  (3002^)  +  1.7471  (300Tq) 

250  mb  =  -  469.23  -  0.97097(3002^)  -  0.90401  (300T^)  +  0.10183  (250T_ ^2) 
+  0.97097  (2502jj) 

200  mb  2q  =  676.82  +  0.96655  (2502^)  +  1.7691  (250Tp) 

200  mb  =  -  317.69  -  0.52450  (2502^)  -  0.33968  (250Tq)  +  0,52450(2002^) 
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150  mb  Zjj  «  718.38  +  1.0007  (2002^^)  +  2.0654  (200Tp) 

150  mb  Tq  =  -  395.51  -  0.51998  (2002^)  -  0.49043  (200Tq)  +  0.51998  (1502^) 

100  mb  Zq  =  588.05  +  1.0119  (1502^')  +  2.3079  (ISOT^)  +  0.39626  (1002_ ^2^ 

-  0.39626  (1502_^2) 

100  mb  =  -  247.42  -  0.1953258  (1502^^)  +  0.26362  (100T_ ^2)  +  0.20314  (1002 q) 

B.2  Kavy  Equations,  for  50  and  30  mb 
B.2.1  January 
B.2.1.1  0-19*N 


50  mb 

^0 

=  2541.416  +  0.819(1002^)  +  2.967  (lOOT^) 

50  mb 

^0 

=  41.201  -  0.015  (1002^)  +  0.291  (lOOT^) 

30  mb 

=  677.860  +  1.082(502^)  +  2.859  (SOT^) 

30  mb 

^0 

=  -  234.379  +  0.030  (50Zq)  +  0,422  (SOT^) 

B.2.1.2 

20-29"N 

50  mb 

=  3533.167  +  0.619  (IOOZq)  +  1.697  (lOOT^) 

50  mb 

^0 

=  101.667  -  0.030  (1002^)  +  0.002  (lOOT^) 

30  mb 

=^0 

=  411.119  +  1.122  (502 jj)  +  2.919  (50Tq) 

30  mb 

■"o 

=  -  296.825  +  0.040 (502jj)  +  0.418 (SOT^) 

-  B.2.1.3 

30-39‘’N 

50  mb 

=  2197.972  +  0.901  (1002^)  +  4,206  (lOOT^) 

50  mb 

^0 

=  -  16.853  -  0.003  (1002jj)  +  0.421  (lOOT^) 

30  mb 

=  1006.255  +  1.045  (502q)  +  4.210  (SOT^) 

30  mb 

^0 

=  -  111.625  +  0.015  (502jj)  +  0.785  (SOT^) 

B.2. 1.4 

40-49'N 

50  mb 

=  1357.186  +  1.074  (1002q)  +  5.453  (lOOT^) 

50  mb 

^0 

=  -  161.097  +  0.027  (1002q)  +  0.722  (lOOT^) 

30  mb 

=  1156.571  +  1.028  (502q)  +  4.948  (507^) 

30  mb 

^0 

=  -  37.062  +  0.006  (502q)  +  1.060  (50Tq) 
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B.2.1.5  SO-SyN 


50  mb 

=^0 

=  1703.937  +  1.027(1002^)  +  7.247  (lOOT^) 

50  mb 

^0 

=  -  39.900  +  0.009  (IOOZq)  +  1.182  (lOOT^) 

30  mb 

"o 

=  1373.147  +  0.997  (502^)  +  5.189  (50Tq) 

30  mb 

^0 

=  24.697  -  0.003(502^)  +  1.104 (50Tq) 

B. 

2.1.6 

60-69”N 

50  mb 

"o 

=  1849.552  +  1.001(1002^)  +  7.448  (lOOT^) 

50  mb 

^0 

=  25.856  -  0.003(1002^)  +  1.230  (lOOT^) 

30  mb 

"o 

=  1399.180  +  0.992  (502^)  +  5.040  (50Tq) 

30  mb 

^0 

=  67.613  -  0.009  (50Zq)  +  1.142(50Tq) 

B. 

2.1.7 

70-90'N 

50  mb 

"o 

=  1776.246  +  1.012(1002^)  +  7.208  (lOOT^) 

50  mb 

^0 

=  5.698  +  0.000  (IOOZq)  +  1.169  (lOOT^^) 

30  mb 

"o 

*  1761.325  +  0.943(502^)  +  5.770  (501^) 

30  mb 

^0 

=  177.576  -  0.024  (50Zq)  +  1.326  (SOT^) 

B.2.2 

April 

B.2.2.1 

O-WN 

50  mb 

^0 

=  2884.129  +  0.761  (1002^)  +  3.233  (lOOT^) 

50  mb 

^0 

=  213.412  -  0.046(1002^)  +  0.372  (lOOT^) 

30  mb 

"o 

=  514.297  +  1.103(502^)  +  2.495  (50Tp) 

30  mb 

^0 

=  -  223.826  +  0.027(502^)  +  0.253  (50T^) 

B.2.2.2 

20-29”N 

50  mb 

"o 

=  2757.580  +  0.785  (1002^)  +  3.214  (lOOT^) 

50  mb 

^0 

=  63.839  -  0.021  (1002^)  +  0.136  (lOOT^) 

30  mb 

"o 

=  409.690  +  1.124  (502q)  +  3.069  (50T^) 

30  mb 

^0 

=  -  271.813  +  0.037  (502^)  +  0.473  (50T^) 
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B.2.2.3  30-39'N 

50  mb  Zq  =  2253.952  +  0.891  (lOOZ^^)  +  4.094  (lOOT^^) 

50  mb  Tq  =  13.654  -  0.009  (lOOZ^)  +  0.355  (lOOT^) 

30  mb  Zq  =  768.846  +  1.079  (50Zjj)  +  4.040  (50T^) 

30  mb  Tq  =  -  217.029  +  0.030  (SOZ^^)  +  0.728  (50Tq) 

B.2.2.4  40-49*N 

50  mb  Zq  =  1992.518  +  0.940  (lOOZ^^)  +  3.999  (lOOT^^) 

50  mb  Tq  =  -  14.231  -  0.003  (lOOZ^)  +  0.429  (lOOT^) 

30  mb  Zq  =  999.604  +  1.046  (50Zq)  +  4.255  (50T^) 

30  mb  =  -  111.154  +  0.015  (50Z^)  +  0.814 (50Tp) 

B.2.2.5  50-59*N 

50  mb  Zq  =  1738.520  +  0.999  (lOOZ^)  +  4.978  (lOOTjj) 

50  mb  Tq  =  -  35.693  +  0.033  (lOOZ^)  +  0.622  (lOOT^) 

30  mb  Zq  =  1093.054  +  1.032  (50Zjj)  +  4.280  (SOT^^) 

30  mb  Tq  =  -  91.110  +  0.012  (50Z^j)  +  0.828  (50T^) 

B.2.2.6  60->69°N 

50  mb  Zq  =  1565.265  +  1.045  (lOOZ^)  +  6.320  (lOOTjj) 

50  mb  =  -  69.546  +  0.012  (lOOZ^)  +  0.903  (lOOT^j) 

30  mb  Zq  =  1011.577  +  1.047  (50Zjj)  +  4.681  (50T^) 

30  mb  Tq  =  -  128.145  +  0.018(50Zjj)  +  0.915(50T^j) 

B.2.2.7  70-90'’N 

50  mb  Z^  =  1521.312  +  1.060  (lOOZ^)  +  7.084  (lOOT^) 

50  mb  =  -  73.812  +  0.015  (lOOZ^)  +  1.160  (lOOT^) 

30  mb  Zp  =  1151.790  +  1.028  (50Z^)  +  4.947  (SOT^) 

30  mb  +  -  60.790  +  0.009  (50Z^)  +  1.040  (50T^) 

B.2.3  July 
B.2.3.1  O-IO’N 

50  mb  Zq  =  1329.789  +  1.038  (lOOZ^)  +  2.268  (lOOT^) 
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50  mb  Tq  “  -  130.070  +  0.015  (lOOZ^^) 
30  mb  Zq  =  1353.682  +  0.979  (50Zq) 

30  mb  Tjj  =  -  245.159  +  0.030  (50Zjj)  h 

B.2.3.2  20-29'’N 

50  mb  Zq  =  1888.503  +  0.944  (lOOZ^] 

50  mb  =  -  79.342  +  0.006  (lOOZ^j)  H 

30  mb  =  907.650  +  1.047  (OOZ^)  H 

30  mb  Tq  =  -  111.434  +  0.012  (SOZ^)  h 

B.2.3.3  30-39°N 

50  mb  Z^  =  1688.097  +  0.993  (lOOZ^^: 

50  mb  Tq  =  -  87.255  +  0.009  (lOOZ^)  n 

30  mb  Z^  =  798.473  +  1.069  (50Z^)  h 

30  mb  Tjj  =  -  147.468  +  0.018  (SOZ^)  - 

B.2.3.4  40-49''N 

50  mb  Zq  =  1948.918  +  0.947  (lOOZ^ 

50  mb  Tq  (IOOZq)  - 

30  mb  Z^  =  875.143  +  1.058  (SOZ^)  - 

30  mb  Tjj  =  -  167.540  +  0.021  (SOZ^^)  ■ 

B.2.3.5  50-59'’N 

50  mb  Zq  =  1550.912  +  1.037  (lOOZp 

50  mb  Tq  =  -  101.673  +  0.015  (lOOZ^) 

30  mb  Zq  =  781.138  +  1.073  (OOZ^)  • 

30  mb  Tq  =  -  182.025  +  0.024  (SOZ^)  • 


B.2.3.6  60-69°N 


50  mb 

o 

II 

1454.129  +  1.060  (IOOZq 

50  mb 

^o  =  - 

108.563  +  0.018  (IOOZq) 

30  mb 

^0  = 

1040.348  +  1.042  (50Zq) 

30  mb 

1 

11 

O 

H 

84.710  +  0.012  (50Z^)  + 

+  0.187  (lOOTjj) 

+  2.446  (50Tq) 

K  0.261  (50Tq) 

)  +  2.796  (IOOTq) 
h  0.181  (lOOT^) 
h  2.703  (50Tq) 

I-  0.401  (50T^j) 

)  +  3.597 (IOOTq) 
I-  0.281  (lOOT^) 
h  3.406  (50Tq) 
h  0.499  (50T^) 

)  +  3.613 (IOOTq) 
h  0.324  (lOOT^) 
f  3.421  (50Tq) 

^  0.516  (50Tp) 

)  +  5.176  (lOOTjj) 
+  0.575  (IOOTq) 
f  3.607  (50T^j) 

+  0.657  (50Tq) 

)  +  5.664  (IOOTq) 
+  0.756  (IOOTq) 

+  4.590  (50T^) 
0.950  (50Tq) 
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B.2.3.7  70-90*N 


Zq  =  1691.480  +  1.011  (lOOZ^)  +  4.947  (lOOT^) 

Tq  =  -  17.860  +  0.000  (IOOZq)  +  0.551  (lOOT^) 
Zq  =  1109.365  +  1.026  (50Zjj)  +  3.622  (507^^) 

'Tq  “  ■  35.241  +  0.003  (50Zq)  +  0.633  (507^) 

B.2.4  October 
B.2.4.1  0-19*N 

50  mb  =  2471.473  +  0.845  (lOOZ^)  +  3.594(1007^) 

50  mb  =  -  4.789  -  0.006  (lOOZ^)  +  0.315(1007^) 

30  mb  Zq  =  414.772  +  1.119(502)  +  2.678(507^) 

30  mb  Tq  =  -  256.894  +  0.034(50Z^)  +  0.402(507^) 

B.2.4.2  20-29‘’N 

50  mb  Zq  =  2172.734  +  0.895(100Z^)  +  3.155(1007^) 

50  mb  =  -  28.805  -  0.003  (lOOZ^)  +  0.205  (1007^) 

30  mb  Zq  =  840.373  +  1.056  (50Zjj)  +  2.621(507^) 

30  mb  Tp  =  -  134.085  +  0.015  (SOZ^)  +  0.375(507^) 

B.2.4.3  30-39'’N 

50  mb  Zq  =  1908.796  +  0.949  (lOOZ^)  +  3.501  (1007^) 

50  mb  Tg  =  -  74.893  +  0.006  (lOOZ^)  +  0.251(1007^) 

30  mb  Z^  =  974.336  +  1.044  (50Z^)  +  3.530  (507^) 

30  mb  =  -  103.101  +  0.012  (50Z^)  +  0.565  (507^) 

B.2.4.4  40-49*N 

50  mb  Zq  =  1564.083  +  1.024  (lOOZ^^)  +  4.398(1007^) 

50  mb  Tq  *  -  111.220  +  0.015  (lOOZ^)  +  0.461  (1007^) 

30  mb  Z^  =  861.734  +  1.064  (OOZ^)  +  4.006(507^) 

30  mb  Tq  ®  -  158.949  +  0.021  (50Z^)  +  0.704(507^) 


50  mb 
50  mb 
30  mb 
30  mb 
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B.2.4.5  50-59'N 

50  mb  Zjj  =  1674.144  +  1.016  (lOOZ^)  +  5.560  (lOOT^j) 

50  mb  Tq  =  -  6X.565  +  0.009  (lOOZ^j)  +  0.772  (lOOT^j) 

30  mb  Zq  =  1091.373  +  1.034  (50Zq)  +  4.555 (50Tq) 

30  mb  Tjj  =  -  88.217  +  0.012  (OOZ^)  +  0.898  (50Tq) 

B.2.4.6  60-69'’N 

50  mb  Zq  =  1535.147  +  1.050  (lOOZjj)  +  6.341  (IOOTq) 

50  mb  Tq  =  -  67.682  +  0.012  (lOOZ^)  +  0.971  (lOOT^) 

30  mb  Zq  =  1233.344  +  1.017  (50Zjj)  +  5.103  (50T^j) 

30  mb  Tq  =  -  38.939  +  0.006  (50Zp)  +  1.067  (SOT^) 

B.2.4.7  70-90'’N 

50  mb  Zq  =  1461.713  +  1.060  (lOOZp)  +  6,030  (lOOTp) 

50  mb  Tq  =  -  88.606  +  0.015  (lOOZ^)  +  0.901  (lOOT^) 

30  mb  Zq  =  1319.242  +  1.005  (50Zq)  +  5.232  (50Tjj) 

30  mb  Tq  =  27.912  -  0.003(50Z^j)  +  1.181  (50Tq) 
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APPENDIX  C 

RESULTS  OBTAINED  FROM  VERTICAL-EXTRAPOLATION  EQUATIONS, 
FOR  DEPENDENT  AND  INDEPENDENT  DATA 


•H,  T,  and  Z  ara  thicknass  of  tayar  bataaan  constant-prassura  surfacas,  tamparatura  at  constant. prassura  surfaca,  and  iiaight  of  constant- 
prassura  surfaca,  raspactivaly.  Subscript  0  dasignatas  obsarvation  tima,  subscript  -12  designates  12  hr  before  observation,  and  subscript  n 
designates  a  naan  over  a  layer. 

t Fifty  independent  cases  for  each  pradictand. 
tTvanty-fiva  cases  for  each  pradictand. 


T,  and  Z  arm  thickness  of  layer  between  constant-pressure  surfaces,  tatnperature  at  constant-pressure  surface,  and  height  of  constant- 
pressure  surface,  respectively.  Subscript  0  designates  observation  time,  subscript  -12  designates  12  hr  before  observation,  and  subscript  m 
designates  a  naan  over  a  layer. 

t Fifty  Independent  cases  for  each  predictand. 

{Twenty-five  cases  for  each  predictand. 


*H,  T,  and  Z  ara  thickness  of  layer  between  constant-pressure  surfaces^  temperature  at  constant-pressure  surface,  and  height  of  constant* 
pressure  surface^  respectively.  Subscript  0  designates  observation  tin«,  subscript  -12  designates  12  hr  before  observation,  and  subscript  ' 
designates  a  mean  over  a  layer. 

tFifty  independent  cases  for  each  predictand. 

ITwenty-five  cases  for  each  predictand. 


i 

I 


I . 


h 

n 

X 

I 


n 

o  -o 
o 


■D  ^ 
C  « 
«  > 


3b 


C'« 

ft 

M  C 


*0  -o 


m  9  o  ■o 

r  k 

o  •  •  L 
>H  W  0  £ 

?51sSi 

“  j( "  8l 

X  a  k 

S'-  “Is 

■5f|; 

Sb  •■S'T 

,*8>.^ 

fS  AH.  K 

f 


} 


I 


\ 


\ 

1 


I 

i 


115 


•■H,  T,  snd  Z  ar#  thickness  of  l«?y*r  botwoon  constant-pressur#  surfacas^  tamparatura  at  constant-prassure  surfaca,  and  halght  of  constant- 
prassura  surfaca,  raspacti valy.  Si-bscript  0  dasignatas  observation  tima,  subscript  -12  designates  12  hr  before  observation,  and  subscript  -n 
designates  a  mean  over  a  layer. 

t Fifty  independent  cases  for  each  pradictand. 

|Twantv-fiva  cases  for  each  pradictand. 
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•»H,  T,  and  2  ara  thickness  of  layer  between  constant- pressure  surfaces,  temperature  at  constant-pressure  surface,  and  height  of  constont- 
pressu-e  surface,  respectively.  Subscript  0  designates  observation  time,  subscript  -IS  designates  12  hr  before  observation,  and  subscript  m 
designates  a  mean  over  a  layer. 

t Fifty  independent  cases  fer  each  predictand. 
tTwenty-five  cases  for  each  predictano. 


T,  and  Z  ara  thicknasa  of  layar  batwaan  constant- praaaura  surfacas,  taniparatura  at  constant-prassura  surtaea,  and  haight  of  conttan 
prassura  surfaca,  raspactivaly.  Subscript  0  dasignatas  obsarvation  tina,  subscript  -^Z  designates  12  hr  bafora  obsarvation,  and  subscript 
designates  a  aaan  ever  a  layar. 

t  Fifty  indapandant  cases  for  each  ^adictand. 
tTuanty-fiva  cases  for  each  pradictand. 


